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Preface 


The third High-Speed Research Sonic Boom Workshop was held at NASA 
Langley Research Center on June 1-3, 1994. The purpose of this workshop was to 
provide a forum for government, industry, and university participants to present and 
discuss progress in their research. The workshop was organized into sessions dealing 
with atmospheric propagation; acceptability studies; and configuration design, analysis, 
and testing. Attendance at the workshop was by invitation only. 

The workshop proceedings include papers on atmospheric propagation and 
acceptability studies. Papers discussing atmospheric effects on the sonic boom waveform 
addressed several issues. It has long been assumed that the effects of molecular 
relaxation are adequately accounted for by assuming that a steady-state balance between 
absorption and nonlinear wave steepening exists. It was shown that the unsteadiness 
induced by the nonuniform atmosphere precludes attaining this steady state. Further, it 
was shown that the random atmosphere acts as a filter, effectively filtering out high 
frequency components of the distorted waveform. Several different propagation models 
were compared, and an analysis of the sonic boom at the edge of the primary carpet 
established that the levels there are bounded. Finally, a discussion of the levels of the 
sonic boom below the sea surface was presented. Papers on acceptability examined the 
reaction of people to long-term sonic boom exposure. 


The use of trade names of manufacturers in this report does not constitute an official endorsement 
of such products or manufacturers, either expressed or implied, by the National Aeronautics and 
Space Administration. 


David A. McCurdy 
Langley Research Center 
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AN OPTICAL TECHNIQUE FOR EXAMINING AIRCRAFT 
SHOCK WAVE STRUCTURES IN FLIGHT 


Leonard M. Weinstein 
NASA Langley Research Center 
Hampton, VA 
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INTRODUCTION 


The detailed properties of sonic booms have to be better understood before commercial, next 
generation, supersonic and hypersonic aircraft can be properly developed. Experimental tests and 
measurements are needed to help sort the physical details of the flows at realistic test conditions. 
Some of these tests can be made in wind tunnels, but the need for full flight conditions simulation, 
the problem of tunnel wall interference, and the short distance the shocks can be examined from the 
aircraft, limit the usefulness of wind tunnel tests. 

Previous measurement techniques for examining the flow field of aircraft in flight have included 
pressure measurements on the aircraft, ground based pressure measurements, and flow field 
measurements made with chase aircraft. Obtaining data with chase planes is a slow and difficult 
process, and is limited in how close it can be obtained to the test aircraft. A need clearly existed for 
a better technique to examine the shock structure from the plane to large distances from the plane. 

A new technique has been recently developed to obtain schlieren photographs of aircraft in 
flight (SAF). Preliminary results have been obtained, and the technique holds promise as a tool to 
study the shape and approximate strength of the shock wave structure around the test aircraft, and 
examine shock wave details all the way from the aircraft to near the ground. The current paper 
describes this approach, and gives some preliminary test results. 
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Figure 1 . Focusing schlieren of shuttle model in 0.3 m TCT. 


WIND TUNNEL SCHLIEREN PHOTOGRAPH 


The photograph in figure 1 was obtained with a focusing schlieren system in the 0.3 m 
Transonic Cryogenic Tunnel. It illustrates a failed attempt to obtain a Mach 1.2 flow field over a 

large shuttle model. The interaction between the body shock and the tunnel wall resulted in a 
partially started flow, with both supersonic and subsonic regions. This picture demonstrates the 
difficulty of testing reasonably large models in supersonic wind tunnels. The large amount of flow 
detail visible also shows the value of direct flow visualization. 
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Figure 2 shows an F-l 1 1 flying just over Mach 1, with a non uniform background sky. The 
strong shock waves distort the image of the background sky by refraction, enough to displace 
features and make the shocks visible. The low sensitivity of this technique, and the need for the 
non uniform background greatly limits this direct imaging technique. 
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Figure 2. Refraction image of F-l 1 1 flow field 
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FOCUSING SCHLIEREN 


The new flow visualization technique, which will be described later, uses some of the features 
of focusing schlieren systems. A brief description of focusing schlieren is given to help understand 
the new concept. 

The features of a focusing schlieren system are shown in figure 3. An array of points or lines 
from a source grid are imaged by a camera lens. A photographic negative image of the source grid 
is then obtained to use as a cutoff grid. Each point on the source grid now has a corresponding 
"knife edge" spatial filter, and each source and cutoff point pair make an uncollimated conventional 
schlieren system with a very small field of view. However, the multiple sources from the source 
grid are selected to be close enough together so that each point on the final test section image is 
obtained from several source points. The net result is that a large field of view is obtained, with 
different parts of the source grid being used to obtain different parts of the image. The field of view 
is now determined by the total size of the source grid (rather than lens diameter). 


Fresnel Camera lens Image 



Figure 3. Details of Focusing schlieren system. 
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SCHLIEREN FOR AIRCRAFT IN FLIGHT SETUP 


The schlieren for aircraft in flight (S AF) was developed by combining some of the features of 
astronomical photography with focusing schlieren and streak photography. The setup used is 
shown in figure 4. A telescope is used to track the sun (or moon) so that the image is stationary at 
the focus. A neutral filter is used to lower intensity to a safe level. An opaque mask with a narrow 
curved slit is positioned to block all but a thin portion of the edge of the sun and a small region of 
sky next to the edge. A thin region at the edge of the sun is thus the light source, and the mask is 
the cutoff, for a single curved grid line as in a focusing schlieren system. An aircraft is flow 
through a field of view crossing this sliver of light, and the aircraft image is sharply focused on to a 
film plane. The film is moved to follow the aircraft image movement. A narrow slice of a schlieren 
image of the aircraft flow field is formed at each instant of time, and over the full exposure this slice 
is scanned along the image. This produces a composite streak camera schlieren image of the flow 
field examined. 



Figure 4. Setup of schlieren for aircraft in flight. 











RESOLUTION AND FIELD OF VIEW FOR SAF 


The maximum possible resolution of detail in images is affected by several factors. The 
maximum resolution (diffraction limit) of a telescope is determined by the telescope aperture. In 
addition, the atmospheric turbulence, accuracy of focus, tracking accuracy during exposure, and 
film limitations may result in less than diffraction limited resolution. Unless the telescope is on a 
high mountain, or in a chase plane, the daytime atmospheric limited resolution is seldom better than 
1 arcsec, and more often closer to 2 arcsec. To prevent significant additional limitation due to the 
telescope, a telescope of at least 5 inch aperture is desired (which can approach 1 arcsec resolution 
with near diffraction limited performance). For the following table, a net image resolution of 2 
arcsec was assumed for all causes combined. 

The field of view photographed is the physical height at the aircraft distance corresponding to 
the angle of the portion of the sun used for a grid line. The sun covers slightly over 1/2 degree, and 
about 1/3 degree is used for the grid line in SAF. This angle was used to compute the field of view 
in the table. 


Table I. Resolution and Field of View as a Function of Slant Range 


SLANT RANGE 
(ft.) 

(for 2 arcsec) 
RESOLUTION 
(inches) 

( for 1/3 deg) 
FIELD OF VIEW 
(ft.) 

8,000 

0.9 

46 

16,000 

1.9 

92 

32,000 

3.7 

184 

64,000 

7.4 

368 

128,000 

14.9 

736 
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PHOTOGRAPH OF SETUP AT WALLOPS FLIGHT FACILITY 


The prototype system used to demonstrate the SAF concept used an 8 inch aperture, f/10 

telescope, which was large enough to carry a reasonable amount of equipment, and was still small 

enough to carry around and set up easily. 

The SAF system was carried to Wallops Flight Facility and set on the end of a runway. Figure 
5 shows the system set up for the initial effort. Several tests flights were made with a T-38 at 
M=0.9 (but with supersonic flow over the wings) at a slant range of 16,000 ft., to try and 
photograph the flow field at high resolution. The flights were required to be subsonic for these tests 
since they were made over land. The field of view was only 46 ft. at that range. A combination of 
clouds, equipment problems, and the difficulty of passing through this small area prevented 
success in the early flights. 



Figure 5. SAF set up on runway at Wallops. 
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SAF OF T-38 AT 32,000 FT. AND M=l.l 


The aircraft range was increased to 32,000 ft. and the plane was flown fully supersonic for 
additional tests. This resulted in a 92 ft. field of view, and a large extent of supersonic flow. 

Assateague island (close to Wallops) was selected as the site to set up the SAF equipment. This 
location was used because the supersonic flights had to be made over water, 5 miles away from 

land, to prevent sonic booms from bothering people. 

At 10:30 a.m.on Dec. 13, 1993, the T-38 was flown at M=l.l at an altitude of 13,700 ft. This 
resulted in a 32,000 ft. slant range from the camera. The Wallops ILS system was used to guide the 
plane through the test area. The photo shown in figure 6 was the first SAF photo of the T-38 flow 
field. 

The image is very grainy because of the film used for the preliminary test. In addition, the relay 
lens used to reimage the aircraft onto the film was not well corrected except near the center of the 
image. Since the plane was near the edge of the field of view, the image was noticeably less sharp 
around the plane. In addition, some banding (curved streaks due to uneven film motion) resulted in 
uneven exposure of the image. Even with these limitations, the image clearly showed the shock 
wave details. About one plane length above the aircraft, 6 shocks are seen to merge into 4. 



Figure 6. SAF of T-38 at 32,000 ft. slant range and M=l.l. 


USE OF SAF TO STUDY SONIC BOOMS 


The demonstration of SAF as a tool to examine the flow field of aircraft in flight led to the 
possibility of using it as part of a sonic boom study. A joint cooperative effort was initiated at 
Langley to evaluate the possible use of SAF to examine both near and far field shock structure for 
aircraft used in these studies. 

A shock wave system weakens as it moves away from a supersonic aircraft, and the shock 
system also distorts as it moves through local variations in the atmosphere. In addition, the aircraft 
used in sonic boom studies are generally flown at significantly higher altitudes than the T-38 flight 
was. This means that lower atmospheric densities are encountered, and the resolution is reduced for 
more distant subjects. The SAF system was taken out to Dryden Flight Facility to evaluate these 
effects, and see if limitation would be encountered. 

Modifications were made to overcome some of the problems encountered in the Wallops test. A 
finer grain film was selected to cut graininess. The relay optical system was modified to obtain 
better sharpness on the outer part of the image. The film used was also selected to have higher 
contrast, to increase sensitivity. This last feature also made the "banding" much worse, but it was 
thought this would not prevent evaluation of the performance, even though it would make bad 
looking pictures. 

The need to fly at fairly large altitudes combined with the desire to keep slant range reasonably 
small, and the need for the atmosphere to be as undisturbed as possible, constrained the flights to 
mid morning. 
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The SAF system was set up on the peak of Galileo Hill in California City. This is close 
to Dryden, and was in a good location to obtain SAF pictures. Figure 8 shows the hill, 
which is about 700 ft taller than the surrounding area. 


Figure 8. Photograph of Galileo hill. 
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Figure 9. SAF of F-18 at 60,000 ft. slant range and M=1.4 


SAF OF F-18 AT 60,000 FT RANGE AND M=1.4 


On April 13,1994, the first SAF test was made at Dryden. Two F-18's were flown at 
M=1.4, one mile apart, with the trailing plane 300 ft. below the leading plane. This was 
done to increase the chance of photographing at least one of the planes. The leading plane 
was about 35,000 ft. above sea level, and the SAF system was about 3,200 ft. above sea 
level. This resulted in a slant range of about 60,000 ft. from the planes to the camera at 
closest approach. The viewing angle of the sun results in an oblique view rather than a side 
view of the aircraft. The photograph shown in figure 9 was taken just after 9 a.m., and 
shows the flow field around the lead plane. 

SAF is sensitive to density gradients normal to the local grid line. Since the grid line 
curves, and since the shocks swept back, the best orientation of grid is to curve back 
toward the shock direction. The photo in the figure had the grid reversed, and this resulted 
in lower sensitivity, especially near the edges of the photo. 

The camera was also slightly out of focus, so the resolution was lower than expected. 
The image banding was very bad, as was expected, but did not prevent the flow field from 
being seen. The aircraft shape shows up well, and shocks are clearly visible, even though 
the plane was about 12 miles away. 
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,060 ft.below F-18. 


On April 15, several additional flights were made by the two F-18's to try to get shock 
detail at different distances below the aircraft. The grid was reversed from- the previous 
photo to obtain the best sensitivity possible. The photo shown in figure 10 was obtained 
just before 9 a.m., for the trailing plane. The photo shows the shock at 1,060 ft. below the 
plane. The complex multiple shock pattern has already merged to just two visible shocks. 


. BELOW F-18 
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VARIATION OF LONGITUDINAL SHOCK SEPARATION 

WITH DISTANCE 


Three passes by the flights made on April 15 resulted in 6 shock photos being obtained. 
These ranged from 510 ft. below the plane (diagonal distance) to 3,830 ft. below. The 
photographs for the shocks at the largest distance were not correctly exposed, so were not 
as good as the others. However, all shock photos were usable. The maximum distance the 
shocks can be photographed has not been determined, but clearly is very far from the 
aircraft. Figure 1 1 shows the change in length of the separation of front to back shock, 
normalized by the plane length, plotted against the diagonal distance below the aircraft. 
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PROBLEMS FROM PROTOTYPE SAF, AND 
SUGGESTED SOLUTIONS 


There were several problems encountered in the prototype SAF system. A list of these, 
and possible ways to overcome them (including some that have been already implemented), 
are given in table 2. 


Table II. Problems From Prototype SAF System, and Suggested Solutions 


PROBLEMS 


SOLUTIONS 


Poor focus with low power eyepiece 
Grainy and low resolution film 

Film banding 

Sun tracking errors 

Camera rotation with equatorial tracker 
Fly plane through correct location 


Use higher power eyepiece, or auto focus 

Fine grain, high resolution film, or high 
resolution electronic imagery 

Smooth film transport, or TDI electronic 
imaging 

Auto sun tracker 

Use altazimuth rather than equitorial tracker 

Use precision navigation, or multiple SAF 
cameras with overlapping fields of view 



ALTERNATE VERSIONS OF SAF 


The prototype SAF camera was used as a development system, and was far from 
optimized. Different versions can be made to better accomplish specific tasks. Table 3 
shows some alternate versions and what special capabilities can be accomplished with 
them. 


Table in. Variations of SAF Systems 
VARIATION SPECIAL CAPABILITY 


Chase plane version 

Small compact system 

TDI camera instead of film 

SAF multiple camera array (one station) 

SAF multiple station array 

Point SAF camera with photo detector 

Near ground source/detector system 


High resolution, possibility of near side view 

Easy portability, flight use 

No banding, instant replay 

Large field of view 

Flight Dynamics 

Time history of flow field at single point 
Near ground time history of signal 



SUMMARY 


The basic capability of schlieren for aircraft in flight imaging technology has been 
demonstrated at NASA's Wallops Flight Facility, and Dryden Flight Facility. Photos of T- 
38 and F-18 aircraft along with their flow fields, have been obtained at 6 and 12 miles 
range respectively, and shock wave pictures have been obtained as far as 3,830 ft. 
diagonally below the F-18. 

Problems in the limitation of the initial design, operational procedures, and recording 
medium have been either solved, or possible solutions suggested. 

An improved SAF system is currently being designed, and will be tested when ready. 
Use of arrays of SAF cameras, and special versions of systems should further increase the 
possible capabilities of this approach. 

The limits of the technology have not yet been reached, but even the demonstrated 
capability should be useful for sonic boom studies. 
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SUMMARY 


Sonic boom propagation is not steady, even in a nonturbulent atmosphere. The shock 
overpressure A p is affected by geometrical spreading, stratification of the atmosphere, and even the 
N shape of the waveform. Nevertheless, for purposes of predicting shock profile and rise time, it has 
commonly been assumed that the shock is in steady state. For example, molecular relaxation, which 
is a major factor controlling sonic boom rise time, is strongly dependent on relative humidity. 
Because humidity varies with altitude, rise time varies as the sonic boom propagates downward. 

The question is whether rise time depends only on local conditions or is also affected by the history 
of humidity variation along the propagation path. Kang [1, Chap. 7.2] argues that shocks respond 
to change in humidity quickly enough that they are in effect always in steady state. In other words 
only local conditions are important. Robinson [2, Chap 5.2] however disagrees with this hypothesis. 
Raspet et. al. [3] found that perturbed 100 Pa shocks (step waveform) require propagation 
distances of order 1 km for the rise time to return to within 10% of its steady shock value. 

f 

The purpose of our investigation is to determine the effect of unsteadiness (not associated with 
turbulence) on rise time. The unsteadiness considered here is due to (1) geometrical spreading, (2) 
stratification, which includes variation in density, temperature, and relative humidity, and (3) N 

1 Work supported by NASA 
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shaped waveform. A very general Burgers equation, which includes all these effects, is the 
propagation model for our study. The equation is solved by a new computational algorithm in 
which all the calculations are done in the time domain. 

The present paper is a progress report in which some of the factors contributing to unsteadiness 
are studied, namely geometrical spreading and variation in relative humidity. The work of Pierce 
and Kang [4], which motivated our study, is first reviewed. We proceed with a discussion of the 
Burgers equation model and the algorithm for solving the equation. Some comparison tests to 
establish the validity of the algorithm are presented. The algorithm is then used to determine the 
distance required for a steady-state shock, on encountering an abrupt change in relative humidity, 
to reach a new steady state based on the new humidity. It is found that the transition distance for 
plane shocks of amplitude 70 Pa is about 4 km when the change in relative humidity is 10%. Shocks 
of amplitude 140 Pa require less distance. The effect of spherical and cylindrical spreading is also 
considered. We demonstrate that a spreading shock wave never reaches steady state and that its 
rise time will be less than the equivalent steady state shock. Finally we show that an N wave has a 
slightly short rise time than a step shock of the same amplitude. 


OUTLINE 

• Introduction 

• Review of Pierce-Kang result 

• Burgers equation model 

• Time-domain algorithm 

• Validation tests of algorithm 

• Rise time as a function of 


1. Change in relative humidity 

2. Geometrical spreading 

3. Waveform 



REVIEW OF THE PIERCE-KANG PREDICTION OF SONIC BOOM RISE TIME 


Pierce and Kang [4] calculated sonic boom rise time by solving a model propagation equation 
called Burgers’ equation [5]. An important contribution was to augment the classical Burgers 
equation with terms that describe molecular relaxation due to both nitrogen and oxygen. By 
assuming that the sonic boom shock near the ground is in steady state, Pierce and Kang simplified 
their Burgers equation and solved it numerically. Rise time predictions obtained from their solution 
were then compared with rise time data from a large number of sonic boom measurements taken at 
Edwards Air Force Base in 1987 [6]. In their initial comparison almost all the measured rise times 
lay well above the predicted curve [7], as shown in the left-hand plot of Fig. 1. The discrepancy was 
attributed to atmospheric turbulence. Later, however, when they made a correction for pressure 
doubling at the ground, they found their prediction to fall in the middle of the data [8]; see the 
right-hand plot in Fig. 1. Although the data would seem to corroborate their prediction, the role of 
turbulence casts doubt on this conclusion. Model experiments done at our laboratory [9, 10] 
indicate that turbulence almost always increases rise time, rarely decreases it. If turbulence does 
have this effect, then molecular relaxation should be expected to provide a lower bound for the 
data, not approximate a mean for the data. Why then does the Pierce-Kang prediction not serve as 
a lower bound? 



Figure 1: The left-hand plot shows Kang and Pierce’s initial prediction [7]. The right-hand plot 
shows the corrected prediction [8]. 
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Our hypothesis to answer this question is based on results found in Refs. [2, 3]. We suspect that 
the shock wave at the head of a sonic boom does not respond quickly enough to variation in 
atmospheric conditions (and to other changes that affect the profile, such as geometrical spreading 
and even wave shape), to justify the steady-state assumption. If our hypothesis is correct, then to 
improve on the Pierce-Kang prediction requires that more than local conditions be taken into 
account. Past history along the propagation path must be significant. Figure 2 shows profiles of 
temperature and relative humidity measured during a sonic boom experiment [1, pp 157-161]. It is 
seen that conditions can change rapidly, particularly during the lower part of the propagation path. 
To determine whether the sonic boom profile can respond quickly to changes of this order, we have 
calculated the effect on rise time of an abrupt change in atmospheric conditions. Thus far we have 
concentrated on changes in relative humidity. For purposes of this paper, temperature and pressure 
are fixed at their ground level values. 




Figure 2: Atmospheric conditions measured during the sonic boom experiment. 
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THE BURGERS EQUATION 


The “classical” Burgers equation is the standard model equation for plane finite-amplitude waves 
in a thermoviscous medium: 


dp f3 dp 2 _ d 2 p 

dx 2p 0 cjj dt' TV dt' 2 ' * 

Here p is acoustic pressure, t time, t' — t — x/cq retarded time, x distance, Co small-signal sound 
speed, po ambient density, /3 coefficient of nonlinearity, and 6tv the thermoviscous loss coefficient 
Pierce [11] added terms to account for relaxation processes. Each relaxation process v is 
characterized by a relaxation time r„ and a change in small-signal sound speed (Ac)j, due to the 
relaxation. In operator notation Pierce’s “augmented Burgers equation” may be written 


d 2 p v (Ac)t/7V fa 

dx 2po4 dt 1 TV dt' 2 V Co l + r »w P 


( 2 ) 


Equation 2 is still for plane waves. If geometrical spreading is included, the equation becomes 


dp a 

di + x P 


(3 dp 2 
2p 0 Co dt' 


— Sty 


d 2 p 

dt ' 2 


+ £ 


U 


(A c)j/jp qp? 
Co 1 + 


( 3 ) 


where the spreading factor o is 0 for plane waves, ^ for cylindrical waves, and 1 for spherical waves. 
This is the equation we have solved numerically to obtain the results reported in this paper. 


Burgers’ equation may be further generalized to include effects of (1) stratification and (2) 
diffraction. Stratification may be included by scaling the dependent variable p and stretching the 
independent range variable x [12]. The Burgers equation for this case is not considered in this 
report but will be taken up in the future. To include diffraction effects, one must use the KZK 
equation, which is a multi-dimensional form of Burgers’ equation; see, for example, Refs. [13, 14]. 
As a spinoff from the present work, relaxation effects have been included in a computer code that 
solves the KZK equation, but no formal report of the results has yet been given. 
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NUMERICAL PROPAGATION 


Solutions of the generalized Burgers equation that are not in steady state involve solving a 
partial differential equation. Except for a few rare cases the solution can only be obtained 
numerically and it is common to use some sort of marching scheme. A time waveform is digitized 
with M samples and then small steps are taken in the propagation direction. At each step 
absorption and nonlinearity are solved in series. It is popular to do the absorption effects in the 
frequency domain as this requires M complex multiplications. However in the frequency domain the 
nonlinear term involves a convolution — which requires of the order of M 2 operations. To speed up 
the code the nonlinear distortion can be applied in the time domain as it requires only order M 
operations. To go between the time and frequency domain one can use the fast fourier transform 
which requires order MlogM operations. Algorithms like the Pestorius [15] code flip-flop between 
the time and frequency domains at each step to take advantage of computing absorption in the 
frequency domain and nonlinear distortion in the time domain, the price being the use of the FFT. 
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TIME DOMAIN ALGORITHM 


It would be nice to stay in one domain but without having to pay the price of M 2 operations for 
a convolution. Lee and Hamilton [13, 14] have developed a method of computing the absorption in 
the time domain. They approximate the absorption with a finite difference equation. This yields a 
tridiagonal matrix system which can be solved in order M operations. The code they developed was 
used to solve the KZK equation, which is a generalization of the Burgers equation to include 
diffraction effects. We have extended the code to account for molecular relaxation and spreading 
effects. The work presented here does not, however, include diffraction effects. 

Apart from its numerical advantage the fully time domain algorithm has the nice property that it 
can propagate pulses. Because the FFT isn’t used, the endpoints of the waveform need not match 
to make a periodic waveform. Therefore step shocks and N waves are easily dealt with. The 
algorithm is particularly suited to finding a steady state solution. Raspet et. al [3] had to use a 
square pulse to find the steady state behavior of a shock. Square pulses have a limited propagation 
range before they turn into sawtooth waves. In the time domain code a pure shock can be 
propagated with out difficulty. 


x x=x+Ax 


r+> 


Thermoviscous Losses 


— ► 

A p(t') = B p(t’) 






Relaxation 


— ► 

A v p(f) = B v p(t') 

-i 





Nonlinearity 


— > 

f = r + -^3 P (r)Ax 
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Figure 4: Time domain approach to solving the Burgers equation. 
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STEADY STATE: NONLINEARITY AND THERMOVISCOUS LOSSES 


A number of cases were run to test the validity of the code. The first was to obtain the 
steady-state solution of the classical Burgers equation for a thermoviscous fluid. The known 
analytical solution for the steady shock is the hyperbolic tangent function. Figure 5 shows how a 
shock front is propagated with the time domain code; a is the distance variable. The first figure 
shows the initial profile, chosen because it looked interesting. The other figures show how the profile 
develops. The final figure, a — 2, shows that the numerical result agrees very well with the 
analytical solution. 
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Figure 5: Propagation of a shock in a thermoviscous medium. 
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STEADY STATE: NONLINEARITY AND ONE RELAXATION PROCESS 


We now determine whether the relaxation part of the code behaves correctly. Polykova et. al. 
[16] obtained the steady-state solution for a finite amplitude wave in a medium with one relaxation 
process but no thermoviscous losses. Their result (denoted PSK in Fig 6) is 


, (l+p/po) D 1 
Tv n (1 +p/Po) D+1 ’ 


where 


D _ (Ac) t ,p 0 co 
Po/3 

Figure 6 shows the result from the propagation program in a monorelaxing fluid. For the values 
chosen relaxation was not enough to stop the waveform from becoming multivalued. In the 
analytical result weak shock theory was used to ensure a single valued function. Multivaluedness 
was prevented in the numerical algorithm by including a small amount of thermoviscous losses. 



Figure 6: The top figure shows the analytical result for the steady-state solution in relaxing medium 
with no thermoviscous effects; D — 0.5. The middle figure shows the initial and steady-state profiles 
obtained by the time domain code. The bottom figure compares the analytical and numerical steady- 
state profiles. 
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STEADY STATE IN AIR 


The last verification test was against Kang’s result. A plane shock front is sent into a standard 
atmosphere with a relative humidity of 10%. This allows us to compare results with a steady-state 
result in Kang’s thesis. In our calculation the shock started out with a hyperbolic tangent profile 
and was then propagated until the profile no longer changed. Figure 7 compares the two results. 



Figure 7: Steady-state solution in air; T — 20 °C, Po = 1 atm, and a relative humidity of 10%. On 
the left is Kang’s profile [1, Fig. 5.8]. On the right is the profile from the time domain code. 
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EFFECTS ON RISE TIME 


Wo will now use the time domain code to investigate the behavior of the rise time of shocks in 
the air. In all calculations the temperature is 20 °C, and the pressure is 1 atm. 

First we examine how long it takes waveform to recover from a small but abrupt change in 
relative humidity. A plane wave shock is propagated in air of given relative humidity until it reaches 
steady state. This steady-state waveform is then used as the input waveform for an atmosphere 
with another relative humidity. 

Second we investigate the effect of spreading on the rise time of a shock front. Shocks that are in 
steady state are propagated as spreading waves; because the amplitude of the shock decreases the 
rise time is expected to increase. We present results for the combined effects of spreading and a 
change in relative humidity. Finally we investigate the difference between the rise time of an N wave 
and a step shock. 


TRANSITION DISTANCE FOR 70 Pa SHOCKS 


We use the term transition distance to describe how far a shock needs to travel to go from one 
steady-state profile to another. A somewhat similar term “healing distance” is commonly used in 
literature related to turbulence for the distance a perturbed shock needs to return to its original 
state [3]. In this case we shall look at transition distances due to a change in relative humidity. 

Figure 8 shows rise time as a function of propagation distance for a plane step shock of 
amplitude 70 Pa which starts in a medium of 20% relative humidity. The relative humidity of the 
second atmosphere is 10%, 20% or 30%. The results show the transition distance to be at least 
5 km. The plot in Fig. 9 shows the rise time for a shock initially in an atmosphere of 50% relative 
humidity. Transition distances are greater than 2 km. 

The initial fluctuations in the rise time are due to rather gross changes in the profile. The 
changes are such that the 10% to 90% definition is not a very suitable measure of rise time. Similar 
fluctuations were observed by Raspet et. al. [3]. 
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Figure 8: Change in rise time for a waveform leaving a medium of 20% relative humidity. 



distance (m) 


Figure 9: Change in rise time for a waveform leaving a medium of 50% relative humidity. 
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TRANSITION DISTANCE FOR 140 Pa SHOCKS 


Figure 10 shows the behavior of a shock wave of amplitude 140 Pa which starts off in air of 20% 
relative humidity. The transition distance needed when the new relative humidity is 30% is about 2 
km. However the transition distance is at least 6 km when the new relative humidity is 10%. 



Figure 10: Change in rise time for a 140 Pa shock leaving a medium of 20% relative humidity. 
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SPREADING 


The amplitude of a spherically spreading shock should decrease as 

Ap = —Ap 0 , (4) 

x 

and for cylindrical spreading, 

Ap = ,/^-Apo • 

V x 

The steady state rise time found from the analytical solution of the classical Burgers equation is 

AT = ln(9) «2gi (5) 

As the amplitude of the waveform decreases the rise time increases because the nonlinear steepening 
effects are weaker. However it is not clear that a spreading waveform will be in steady state. This 
would require the absorption mechanism to respond immediately to the spreading. Naugol’nykh [18] 
argued that a spreading shock in a thermoviscous medium should have a rise time that is shorter 
than the steady-state value because the absorption mechanism can’t work fast enough. 

If a spreading shock remains in steady state then, from Eq. (5) the rise time should vary 
inversely as the pressure jump. Since for spherically spreading waves the pressure varies inversely 
with distance, Eq. (4), we have 

At oc x . 

To investigate the validity of this relation we started with the hyperbolic tangent profile appropriate 
for a plane step shock. The shock was then propagated as a spreading wave. Figure 11 shows the 
initial waveform and how the shock diffuses as it loses amplitude. Figure 12 compares the steady 
state prediction of the rise time to the numerically calculated rise time. Note that the steady state 
prediction always overestimates the rise time. Absorption cannot act quickly enough to diffuse the 
profile before more amplitude decrease, due to spreading, occurs. 
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Figure 11: A shock front that starts off in steady-state in a thermoviscous medium is propagated as 
a spherical wave. The upper plot shows the initial waveform. The lower plot shows the waveform at 
20 times the source radius; the steady-state waveform is also shown. 
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Figure 12: The rise time of a step shock in a thermoviscous; medium. The waveform starts off in 
steady-state and is propagated as either a spherically or cylindrically spreading wave. Two different 
initial amplitudes are used. Rise time is normalized to the initial rise time and distance is normalized 
to the source radius. 
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EFFECT OF SPREADING ON TRANSITION DISTANCE 


We now investigate the combined effects of spreading and a change in relative humidity. A source 
is assumed to be 15 km away and the shock front is propagated at an angle corresponding to a sonic 
cone from an aircraft flying at Mach 2. Cylindrical spreading is used. Figure 13 shows the rise time 
curves and compares them to the plane wave curves taken from Fig. 8. In Fig. 14 we compare the 
rise time curves to those of spreading waves in a homogeneous atmosphere, which are the effective 
asymptotes. We see that the waveform is never in steady state but that the asymptotes are reached 
in 6-8 km. 



Figure 13: The rise time of a shock that is cylindrically spreading; the shock starts from a steady-state 
profile for an atmosphere with 20% relative humidity. The dotted lines show the results for plane 
waves as shown in Fig. 8. Note that the steady-state rise time is not achieved. 



Figure 14: The same curves as in Fig. 13, except that the dotted lines show the results for cylindrically 
spreading waves in a homogeneous atmosphere; these are the asymptotic results. 
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EFFECT OF PROFILE ON THE RISE TIME 


Finally we compare rise time for an N wave with that for a steady-state step shock. Figure 15 
shows the curves for a plane N wave that has propagated a long distance in a medium with a 
relative humidity of 50%. The amplitude of the N wave is 48 Pa. We see that a steady-state shock 
of the same amplitude has a longer rise time. However a steady-state shock of amplitude 50.4 Pa 
has the same rise time. It appears that the profile corresponds to that of a higher amplitude shock 
because the amplitude of the N wave decreases with propagation distance. Just as was seen with 
spreading waves, the absorption cannot keep up with the loss in shock amplitude. 



Figure 15: The solid line is an N wave with duration of 125 ms. It has propagated for 30 km. The 
dashed line is a step-shock of the same amplitude. The dot-dash line is a shock of the same rise time 
as the N wave. This shock has an amplitude of 50.4 Pa. 


36 




CONCLUSION 


Our study indicates that assuming a steady-state profile for the head shock of a sonic boom 
tends to cause an overestimate of the rise time. First, we have found that substantial propagation 
distance is required for a shock to respond to abrupt but small changes in relative humidity. When 
the absorption increases, the actual rise time is always less than the steady-shock value. Second, 
neglect of geometrical spreading results in an overestimate of rise time. Third, modeling the N wave 
as a step shock also raises the rise time, although the effect is small. All of these results are 
consistent with the conclusion that for a decaying shock dissipation is always “struggling to catch 
up with nonlinearity.” 
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INTRODUCTION 


One common method of extrapolating sonic boom waveforms from aircraft to ground is to 
calculate the nonlinear distortion, and then add a rise time to each shock by a simple empirical rule. 
One common rule is the "3 over P" rule which calculates the rise time in milliseconds as three 
divided by the shock amplitude in psf 1 . 

This rule was compared with the results of ZEPHYRUS, a comprehensive algorithm which 
calculates sonic boom propagation and extrapolation with the combined effects of nonlinearity, 
attenuation, dispersion, geometric spreading, and refraction in a stratified atmosphere 2 . It is 
shown here that the simple empirical mle considerably overestimates the rise time estimate. In 
addition, the empirical rule does not account for variations in the rise time due to humidity variation 
or propagation history. 

It is also demonstrated that the rise time is only an approximate indicator of perceived loudness. 
Three waveforms with identical characteristics (shock placement, amplitude, and rise time), but 
with different shock shapes are shown to give different calculated loudness. 

This paper is based in part on work performed at the Applied Research Laboratories, the 
University of Texas at Austin, and supported by NASA Langley. 


RISE TIME VS. PRESSURE AMPLITUDE 


A number of waveform propagation calculations were performed, using the ZEPHYRUS 
algorithm, for the case of an aircraft moving at Mach 1.5 at an altitude of 55,000 feet (16,764 m). 
For simplicity, the initial sonic boom waveform was modeled as an N-wave 200 msec in duration. 
The calculations were performed for a number of initial wave amplitudes and ground-level 
humidities. In each case, the initial waveform was extrapolated to the ground, and the ground- 
level amplitude and rise time of the initial shock was calculated (here, the rise time is calculated as 
10% to 90% of peak amplitude). These results are plotted in Fig. 1. This is compared with the "3 
over P" prediction. 

Three features are immediately evident. First, the "3 over P" method has no humidity 
dependence, which is clearly evident in the ZEPHYRUS results. Second, the "3 over P" method 
considerably overestimates the ZEPHYRUS calculation for every case, with the best agreement for 
very dry air. Finally, the "3 over P" method comes closer to agreement with the ZEPHYRUS 
results for high amplitude waves, with increasing overestimation as the amplitude decreases. 
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RISE TIME AND PROPAGATION HISTORY 


The consistent overestimation is probably due, in part, to the fact that the "3 over P" curve is 
based upon the rise time corresponding to a steady-state step shock, or a well-developed shock in 
which the mechanisms affecting rise time have reached a state of equilibrium for the environment at 
ground level. However, it has been demonstrated by Raspet, Bass, Yao, and Wu 3 , as well as by 
this author 4 , that the typical healing distance of the shocks at ground level is on the order of 
hundreds of meters to kilometers, a longer distance than that for which the ground-level conditions 
have been in effect along the propagation history, and thus the shocks at ground level are not in 
steady-state. If the propagation history is accounted for, as in the ZEPHYRUS result, shorter rise 
times are predicted. 

The angle of incidence of the sonic boom with the ground depends upon the aircraft Mach 
number. Longer rise times, closer to steady-state, would be expected for Mach numbers close to 
one, in which the incident wave is closer to a grazing angle with the ground, and thus the wave is 
at near ground conditions for a longer period in its propagation history. Of course, significant 
refraction effects would be expected in this regime. 


RISE TIME AND SHOCK SHAPE 


It has also been demonstrated previously by this author that, for the same ground-level 
pressure amplitude, the rise time varies considerably as a function of the slope of the waveform 
behind the lead shock 4 . Thus, the shape of the waveform near the shock is also important in 
determining rise time. 

The shape of the shock itself is also of intrinsic importance. Although the shock rise time is 
useful as a rough indicator of the expected loudness, too much reliance on the rise time can be 
misleading. The shape of the shock has a strong effect on the high-frequency content of the sonic 
boom and thus on the loudness. In Fig. 2, three waveforms are shown with identical shock 
characteristics (shock placement, rise time, and amplitude), but with different shock shapes: (1) 
linear step shock, (2) shock interpolated with a hyperbolic tangent, and (3) the ZEPHYRUS 
prediction. The loudness of each waveform was computed using the SIG7CZ code (written by 
Brenda Sullivan at NASA Langley), and the results are displayed in Table 1. 


Table 1. Loudness Results 



dB (PL) 

dBC 

dBA 

Step shock 

93.32 

98.27 

78.91 

Tanh 

92.29 

98.29 

78.70 

ZEPHYRUS 

95.59 

98.29 

82.27 


For the PL and A weightings, the ZEPHYRUS prediction is several dB louder than the other 
predictions. The reason can be seen in Figs. 3 and 4, which display the leading and tailing shocks. 
In the ZEPHYRUS prediction, the dispersion due to the molecular relaxation effects is shown to 
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produce a highly asymmetric shock with a very steep initial phase. In turn, the high frequency 
content corresponding to this steep initial phase produces the increase in loudness. 


CONCLUSIONS 


It has been shown here that simple methods of sonic boom extrapolation may be too crude to 
give accurate predictions, as they ignore the effects of humidity, propagation history, and 
waveform shape. In addition, the commonly used "3 over P" rule considerably overestimates the 
rise time. 

The rise time is a useful indicator of loudness, but may be relied upon too much. It has been 
shown here that the shock shape can affect the loudness of the overall waveform. Simple shock 
interpolations may produce a loudness estimate that is too low by several dB. In the case of the "3 
over P" rule, this underestimation of loudness would be further exacerbated by the rise time 
overestimation. 
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Figure 1. Rise time vs. shock amplitude for 2%, 5%, 10%, 25%, 75%, and 95% 
relative humidity at ground level and the "3 over P" curve 
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Figure 2. Calculated waveforms: (a) linear step shock, (b) hyperbolic tangent, and 
(c) the ZEPHYRUS prediction 
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Figure 4. Trailing shock: (a) linear step shock, (b) hyperbolic tangent, 
and (c) the ZEPHYRUS prediction 
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This is essentially a progress report on a theoretical 
investigation of the propagation of transient waves in a random 
medium. The emphasis in this study is on applications to 
sonic -boom propagation, particularly as regards the effect of 
atmospheric turbulence on the sonic -boom waveform. The analysis 
is general, however, and is applicable to other types of waves 
besides sonic-boom waves. 

The phenomenon of primary concern in this investigation is 
the f inestructure of the wave. The first figure shows what is 
meant by the f inestructure . 


NOTE: Research supported by NASA Langley Research Center 
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1. PRELIMINARIES 


The uppermost sketch shows a typical transient waveform as 
it might look if observed under ideal conditions; i.e. , in the 
absence of scattering, refraction, or any other mechanisms which 
might act to distort the waveform. This is what I call the 
incident wave. 

When a wave propagates in a real medium, such as the 
atmosphere, a random high-frequency structure, called the 
f inestructure, generally apears on the waveform, giving it the 
appearance shown in the bottom sketch. The f inestructure is 
defined formally, at least insofar as the present investigation 
is concerned, as the difference between the observed wave and its 
ensemble average, and generally has an appearance much like that 
shown in the middle sketch. 

The f inestructure is believed to arise as a result of 
scattering and redistribution of wave energy by the random 
inhomogeneities of the medium. The f inestructure is particularly 
important in connection with sonic-boom propagation, since the 
strong high-frequency component that often appears near the front 
of the N-wave may contribute appreciably to the perceived 
noisiness of the sonic boom. 

The f inestructure is most conveniently characterized by its 
statistical properties, such as the variance, the standard 
deviation, the correlation function, and the spectral density 
function. The spectral density function is the quantity that is 
most relevant to the question of annoyance, particularly as 
regards sonic-boom waves, since it contains information on the 
frequency content of the finestructure . The spectral density 
function is the quantity of primary interest in this study. 


< 


Incident Wave 
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2. FORMULATION 


The investigation is based on a relatively simple 
mathematical model . In this model the medium is assumed to be 
one-dimensional, with properties that vary only with the spatial 
coordinate; i.e., the medium is assumed to be independent of 
time. The medium is also assumed to be quiescent, which means 
that only thermal scattering is being considered. (Mechanical 
turbulence, which gives rise to inertial scattering, is ignored.) 
In addition, the medium is assumed to be non-dissipative . 
Non-linear effects are disregarded. 

The starting point of the analysis is the time -dependent 
scalar wave equation in one spatial dimension x. (Note that 
derivatives are denoted by superscripts rather than by 
subscripts.) Here w is the wave function; f is the source term. 

The sound speed c*, which is assumed to be a random function 
of x, is written as shown, where c is a constant reference speed 
and n is a stationary random function of x having zero mean 
and unit variance. (The angle brackets denote an ensemble 
average.) The- paramter e is the standard deviation of the index 
of refraction of the medium, and is assumed to be small . 

It should be pointed out that, in the atmosphere, the 
standard deviation of the acoustic index of refraction is 
typically about one part in one thousand, so that the assumption 
that € is small, which is to say, the assumption that the medium 
is' only weakly inhomogeneous, is appropriate for the study of 
sonic-boom propagation. 

The initial conditions reflect the assumption that the 
medium is at rest prior to the initial disturbance. 


c* 2 w tt — w xx = /; t > 0, — oo < x < + oo 

C*(a0 = r ^-7-r , </u>= 0 , < n 2 >= 1 

1 tjjbyXj 

w(x, 0) = w t ( x, 0) = 0 
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3 . PERTURBATION METHOD 

The assumption of a weakly inhomogeneous medium allows the 
problem to be solved by a perturbation method. To begin, 
assume that the wave function w has an expansion (which may be 
asymptotic) in powers of e . By substituting this 
expansion back into the wave equation, equations for the 
coefficients wo , wj , etc. are derived in the usual way. The 
first two such equations, for wg and wi , are shown. It 
turns out to be sufficient to carry out this procedure only as 
far as the first-order term (the term) . The wave function w 
is then approximated by the first two terms of the expansion, as 
indicated on the last line. 

Instead of specifying the source term f, it turns out to be 

more convenient to specify the zeroth-order wave function; i . e . , 

the Wg term, which corresponds to the incident wave. It is 
written as shown, where h, which is an arbitrary function, 
defines the waveform of the incident wave. 

Writing the incident wave in this form is equivalent 
to specifying a source term that is concentrated at the 
origin, having a time dependence determined by the function h. 

With Wg specified in terms of the function h, the right-hand 
side of the equation for wj is determined, provided that M is ' 

regarded as a known function. That equation can be solved by the 

method of characteristics. That method is well known, and 
so will not be described here. Instead, let's 
look next at how an expression for the correlation 
function of the f inestructure is derived, once .an 
expression for wi has been obtained. 


w(x, t; e) = u/g(m, t ) + ewi(x, t ) + e 2 W 2 (x, t) + ... 


C 2 Wg - wg X = / 

— 2 tt XX n —2 tt 

C Wi — wf = 2c fMWQ 


tug (x,t) = h(t — |x|/c) 

w = wq + ew\ + 0(e 2 ) 
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4 . CORRELATION FUNCTION 


To get a general expression for the correlation function of 
the f inestructure, start with the expression for w and average 
it, noting that the average of wj is zero. (The average of wj is 
zero because it is linear in the random function jj, , which has 
average zero.) Note that the average field differs from the 
incident field by terms of order € 2 . 

Next, subtract to get an expression for w, the fluctuating 
field, defined as the difference of w and its average. The 
quantity w corresponds to the f inestructure. An expression for 
the temporal correlation function of the f inestructure at any 
point x is then obtained by forming, at x, the product 
of w with itself at the two time values t.% and t2 • 

The averaging procedure used here is not the usual ensemble 
averaging. Instead, it's a travel-time-corrected averaging 
procedure that I call asynchronous ensemble averaging. The 
essential idea of the method is that, instead of measuring time 
with respect to some universal reference time, such as the time 
that the wave is emitted by the source, it is measured, for each 
wave in the ensemble, relative to the time that that wave arrives 
at the observer. The advantage of using a travel-time-corrected 
averaging procedure of this type is that it avoids certain 
spurious effects that arise as a consequence of averaging over an 
ensemble of waves that have become dispersed due to variations in 
travel time among the different members of the ensemble. 

A travel-time corrected averaging procedure similar to the 
one used here was used by Allan Pierce in one of his papers on 
sonic-boom propagation. 

I should mention also that asynchronous averaging has the 
effect of renormalizing the perturbation series for the wave 
function, which results in the elimination of some secular terms 
from the expression for the first-order field. 

In order to apply asynchronous averaging in deriving an 
expression for the correlation function of the f inestructure, one 
replaces ti by r + r, where r is the travel time from the source 
(the origin)' to the observation point x, and r is the time 
relative to the arrival time at which the observation is being 
made. The variable t2 is replaced by r +■ r +- s, where s denotes 
the separation time between the two observation times . The 
result is an expression for the correlation function of the 
finestructure, which is denoted by K. Note that terms of order e 6 
have been dropped from the expression for K. 

The spectral density function is then obtained by taking the 
cosine transform of K with respect to s. 


(See next page.) 
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CORRELATION FUNCTION 


w =wq + ew\ + 0(e‘ 2 ) 

< w > — wo 4- 0(e 2 ) 
w— < w >= w = ewi 4- 0(e 2 ) 

< w(x,ti)w(x,t 2 ) > = e 2 < wi(x,ti)wi(x,t 2 ) > +0(e 3 ) 

ti—r + r; t 2 — r + r + s; r > 0, s > 0 

K(x, r, s ) =< w(x, r + r)w(x, r + r + s) > 

= e 2 < wi(x,r + r)w\(x,T +r + s) > 
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5. GENERAL RESULTS 


The result of the analysis sketched above is the general 
expression shown for the correlation function of the 
f inestructure . The function 7 that appears in this formula is a 
temporal correlation function that's related to the spatial 
correlation function of the medium by the change of variable 
indicated. The function h determines the waveform of the 
incident wave, as discussed previously; the primes denote 
derivatives. The lower limit of the second inner integral; i.e., 
f, is defined as shown. 

This formula is valid provided that the propagation path 
length x is much greater than the integral scale, or outer scale, 
L, of the medium. (This condition is satisfied in virtually all 
cases of practical interest . ) Note also that this expression is 
valid only for r positive. When r is negative, i.e., before the 
wave has arrived, the correlation function is zero. 

This formula, which expresses the correlation function of 
the finestructure in terms of the incident waveform and the 
correlation function of the medium (the last two quantities both 
being arbitrary) , is the main result of this investigation. 

The spectral density function of the finestructure is 
obtained by taking the cosine transform of K with respect to s 
(the separation variable) . 

One point about this expression that's worth noting is that 
it is independent of the propagation range x. The reason for 
this is that, according to the theory, by the time the incident 
wave has propagated a distance of the order of an outer scale 
length into the medium the finestructure has become, in a 
statistical sense, fully developed. Beyond this range, i.e., for 
x » L, there is, for all practical purposes, no further 
evolution of the statistical properties of the finestructure. 

The correlation function K is, however, strongly dependent 
on the magnitude of the randomness of the medium, as is shown by 
the presence of the term e 2 on the right-hand side of this 
expression. 

Note also that the correlation function K depends on r (the 
elapsed time since onset) , as well as on s, the separation time. 
The spectral density function of the finestructure will therefore 
also depend on r. The correlation and spectral density functions 
of the finestructure are thus time-dependent functions . This is 
a consequence of the fact that the finestructure of a transient 
signal is generally a non-stationary random process. 

The formula shown expresses the function K as the sum of two 
iterated integrals. Of these, the two inner integrals involve 
only the function h. If that function has a sufficiently simple 
form, then the inner integrals can be evaluated analytically, 
after which only a pair of single integrals is left to deal with. 

One case in which this type of simplification is feasible is 
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when the incident wave has the form of a simple ramp function, 
since then the function h is piecewise constant. The relevant 
calculations have been carried out, and the results are shown on 
the next figure. (This type of simplification is also feasible 
when the incident wave has the form of an N-wave, but the 
calculations for that case have not been done.) 


GENERAL RESULTS 


K(r,s) = ^e 2 ^ 7(0 j h' (rj - £)h' (7/ + s)dqd£ 

+ J 7 (€ + s)J- h' (rj - £)h’ {T))dr}d£ 

7(0 = < H s )H s + 0 >» = p Q cs ^ 

£ = max (O,0> x >> L, r > 0, s > 0 
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6. RAMP FUNCTION 


The sketch shows the function h(t) for the case of a 
ramp-function incident wave. The amplitude P and the rise time 6 
are shown on the sketch, as are the two time values r and r+s at 
which the f inestructure is evaluated in order to form the 
two-point correlation function. The results described here are 
for the case in which the points r and r+s are both behind the 
rise phase of the incident wave, as indicated. This is the case 
of most practical importance. 

As has been mentioned, the two inner integrals in the 
expression for K shown on the previous vu-graph can, in this 
case, be evaluated exactly. The result of that calculation is 
the first equation shown. Note that, although K is written 
generally as a function of both r and s, the right-hand side of 
this expression is in fact independent of r, showing that the 
f inestructure is a stationary stochastic process in this case. 

An expression for the spectral density function of the 
f inestructure, denoted by D, is obtained by taking the cosine 
transform of K with respect to s. The result is given by the 
second equation. The function (5 is the transform of 7 • 

Of these two quantities, the spectral density function; 
i.e., D, is, as has already been emphasized, the quantity of 
primary interest in this study. The expression for this quantity 
will be examined more closely in a moment. First, however, 
let ' s take a quick look at what happens to these two formulas in 
the limit as S goes to zero, i.e., as the ramp function becomes a 
step function. 

The limiting forms of K and D, denoted by Ko(s) and Do (to) , 
respectively, are easy to calculate, and are given by the last 
two equations. The function Kq(s) is then the correlation 
function, while Dq (w) is the spectral density function, of the 
f inestructure for the case in which the incident wave is a step 
function. 

It should be pointed out that no analytical problems arise 
in calculating these limits . There is thus no problem with 
singular behavior of the f inestructure in the limit as the ramp- 
function incident wave becomes a step function. 

These last two formulas show that, when the incident wave is 
a step function, the correlation function of the f inestructure is 
proportional to the correlation function of the medium, and, 
similarly, the spectral density function of the f inestructure is 
proportional to the spectral density function of the medium. 

When the incident wave has the form of a ramp function, 
however, with a non-zero rise time, these simple relations no 
longer obtain. This situation, particularly as regards the 
spectral density function, will be examined in more detail in 
the next section. 
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7. RAMP FUNCTION (cont'd) 

In order to examine some of the implications of the result 
obtained above for the ramp-function finestructure spectrum D, 
it's convenient to write the expression for that quantity in the 
form shown on the first line. In this formula the function Dq is 
the finestructure spectrum for the case of zero rise time, i.e., 
for the step-function incident wave, as described above. The 
function A determines how the finestructure spectrum differs from 
Dq when the rise time 6 is non-zero . 

The function A is sketched in the figure. Note that it is 
effectively zero when its argument is greater than about 2 tt. What 
this means is that (referring to the formula for D) the 
finestructure spectrum is effectively zero at all frequencies for 
which 6u is greater than about 27r, or, what is equivalent, when 
f , the frequency in Hz, is greater than about 1/6 . 

This result shows that, the greater the rise time of the 
incident wave, the less high-frequency energy there is in the 
finestructure. A couple of examples of this effect, involving 
parameter values typical of sonic-boom propagation, are shown: 

When the rise time is one millisecond, there is effectively no 
energy in the finestructure spectrum at frequencies above about 1 
kHz. If the rise time is increased to two milliseconds, then the 
lower limit of the frequency range for which the finestructure 
spectrum is devoid of energy drops to 500 Hz. 

Note, however, that, since the function A is equal to unity 
when its argument is equal to zero, the low-frequency portion of 
the finestructure spectrum is relatively unaffected by changes in 
the rise time. 

Since it ' s generally the high-frequency portion of the 
acoustic spectrum that's most annoying to the human ear, what 
these results imply is that, at least for a ramp-function wave in 
one dimension, increasing the rise time of the incident wave will 
tend to reduce the annoyance associated with the finestructure. 
Whether this result holds in three dimensions as well remains to 
be seen, but if it does it has obvious implications as regards 
the idea of shaping the sonic boom in order to reduce its 
annoyance . 

The results just described could, of course, be expressed 
just as well in terms of time scales, rather than in terms of 
frequency. Expressed in those terms, the results imply that a 
ramp- function wave propagating in a random medium will contain 
very little structure having time scales smaller than the rise 
time. This prediction agrees with observations of sonic-boom and 
other types of transient waves , which rarely show any appreciable 
structure having time scales smaller than the rise time. 
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RAMP FUNCTION (cont'd) 


D(r,u) = A(8u>) D 0 (u>) 


The function A(8uS) r and therefore the f inestructure spectrum, is 
effectively zero for 8u> >2 tt, or / > 1/8. 

Sonic boom: 8 — .001 sec; /> 1000 Hz 

8 = .002 sec; /> 500 Hz 
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8. SUMMARY 


The main results of this investigation can be summarized as 
follows. First, the approach described here has been found to 
yield fairly concise, general expressions for the correlation and 
spectral density functions of the finestructure. Second, the 
results for the case of a ramp -function incident wave agree 
with the observation that the important time scales associated 
with the finestructure of sonic-boom and other types of transient 
waves are generally comparable to, or greater than, the rise time. 
Finally, the results (again for the ramp-function incident wave) 
indicate that an increase in the rise time of the incident wave 
is associated with a reduction in the high-frequency content of 
the finestructure. 

This last result, assuming that it holds in the 
three-dimensional case as well, has implications as regards the 
idea of shaping the sonic boom in order to reduce its annoyance. 

It suggests that increasing the rise time of the sonic boom will 
make it quieter, even in the presence of turbulence. 
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INTRODUCTION 


When finite amplitude (or intense) sound, such as a sonic boom, propagates through a turbulent 
atmosphere, the propagation is strongly affected by the turbulence. The interaction between sound 
and turbulence has mostly been studied as a linear phenomenon, i.e., the nonlinear behavior of the 
intense sound has been neglected. It has been shown that turbulence has an effect on the perceived 
loudness of sonic booms, mainly by changing its peak pressure and rise time. Peak pressure and rise 
time are important factors that determine the loudness of the sonic boom when heard outdoors 
[1, 2]. However, the interaction between turbulence and nonlinear effects has mostly not been 
included in propagation studies of sonic booms. It is therefore important to investigate the influence 
of acoustical nonlinearity on the interaction of intense sound with turbulence. 
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OVERVIEW 


As stated in the introduction, the motivation for the research presented here is an investigation of 
the effect of turbulence on the nonlinear propagation of sonic booms and spark-produced N waves. 
In a previous study [2] Lipkens showed that model experiments are successful in simulating the 
sonic boom propagation through a turbulent atmosphere. A nonlinear propagation model for the 
propagation of sound through turbulence is presented here. Results from a numerical experiment 
that simulates the propagation of spark-produced N waves through turbulence are shown. Finally, 
some conclusions of the effect of turbulence on the nonlinear propagation of N waves are given. 


Overview 


• Motivation: study the effect of turbulence on the nonlinear 
propagation of sonic booms (N waves) 


• Model experiment is successful in simulating the sonic boom 
propagation through atmospheric turbulence 


• Nonlinear propagation model 


• Results 


• Conclusion 


Figure 1: Overview of presentation. 
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MODEL EXPERIMENT 


In Fig. 2 the setup of the model experiment is shown. The discharge of a capacitor across the gap 
between the two electrodes creates a spark. The spark-produced N wave propagates across the 
turbulent field of a plane jet. A wide-band condensor microphone picks up the signal. A centrifugal 
fan blows the air into the settling chamber. The air exits through the narrow slit and forms a plane 
jet. A detailed description is presented in Ref. [2]. 



Figure 2: Model experiment setup. 
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In Fig. 3 several waveforms of the model experiment are shown. The upper left graph presents the 
waveform in absence of turbulence. All other waves have propagated through the turbulence. The 
distortion observed in the waveforms is very similar to that in sonic boom signatures, see e.g., [4]. 
Rounded and peaked waveforms are shown, as well as double and triple peaked waveforms. 








Figure 3: Waveform examples of the model experiment. 


I 
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NONLINEAR PROPAGATION MODEL 


We developed a model that simulates the weakly nonlinear acoustic propagation through 
turbulence. The main assumptions are: (1) linear geometric acoustics is sufficient to propagate the 
waves through the turbulent fields and (2) a nonlinear transport equation governs the propagation 
along the rays. We assume that the turbulence is frozen during passage of the wave. Thus, we 
model the turbulent medium as a set of independent realizations of a random field. The turbulent 
field of each realization is constructed as a sum of a finite number of Fourier modes. The model is 
flexible. Two-dimensional (2D) or three-dimensional (3D), scalar (e.g., temperature) or vectorial 
(e.g., velocity) fields can be constructed. The acoustic propagation model then consists of two parts. 
First, linear geometric acoustics is used to trace rays through each individual realization of the 
turbulent field. Then, a Pestorius’ type algorithm [4] is used to solve numerically the nonlinear 
propagation equation along the rays. 


Nonlinear Propagation Model 


• assumptions: - linear geometric acoustics is used to trace rays 
- nonlinear transport equation along the rays 


• turbulence model: - 2D or 3D, temperature or velocity turbulence 

- realizations constructed as a sum of discrete 
Fourier modes 


• ray tracing through each realization of the turbulence 


• nonlinear transport equation: - transformation into Burgers equation 

- solved by a Pestorius type algorithm 

Figure 4: Nonlinear propagation model. 
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TURBULENCE MODEL 


The turbulent field is homogeneous and isotropic with a zero mean value. In Eq. 1 the construction 
of a turbulent velocity field as a sum of N modes is shown [6]. The direction of the turbulent wave 
vector K l of each mode is random and determined by the random angle Oi (in 2D) (Fig. 5). The 
homogeneity of the turbulent field is assured by the random phase shift fa. The amplitude of each 
mode is determined by a predefined energy spectrum. For the results reported here a Gaussian 
longitudinal correlation function was assumed. The associated energy spectrum is shown in Fig. 6. 

N 

Turbulence model tf(5) = ]P u { (ft*) cos (K* ■ x + fa) (1) 

1 

u i (K i )-K i = 0 



Figure 5: Wave vector geometry of a single Fourier velocity mode in 2D. 

Two-dimensional kinetic energy spectrum 

Gaussian longitudinal correlation function f(r) = exp (— r 2 / L 2 ) 

Two-dimensional energy spectrum E(K) = ^-K 3 L A exp (— K 2 L 2 J A ) 

Amplitude of each Fourier mode | fa |= \J~E(K) AK 



K (nr 1 ) 


Figure 6: 2D energy spectrum for a Gaussian correlation function. 
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LINEAR GEOMETRIC ACOUSTICS 


Equations 2-5 present the well-known formulation of acoustic propagation in the geometric 
acoustics limit [7]. The position vector x (Eq. 2) and the nondimensional wave vector p (Eq. 3) at a 
current point on the ray trajectory are completely determined by the value of s and the initial 
position along the incident wavefront. The local geodesic elements (Eqs. 4-5) govern the evolution 
of the wavefront along each ray and permit the calculation of the cross-sectional area of an 
infinitesimal ray tube. The parameters do and do characterize the initial wavefront. The effect of 
the random field is represented by the index N — cq/c of the medium and the Mach number 
M — v/c of the turbulent flow. Once initial conditions are known, a Runge-Kutta fourth-order 
scheme is used to solve the differential system. The step size ds is a function of the maximum wave 
number value of the turbulent energy spectrum, and equals ds — l/2fc max . It is important to note 
that the description of the turbulent field in terms of Fourier modes allows us to obtain analytically 
all the spatial derivatives in the system. Numerical errors are thus reduced. 


Geometrical acoustics: Ray path equations 


I = 
f = 


dR 

ds 

dQ 

ds 


^(Q-V-Qu) 

~ (r • V ' VN - R • (wm) • V- (vm) • Q) 
-j~ 2 {R-VN)(VN-(VM).p) 


( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


where: 

N is the index of the medium, i. e. N = co/c 

M = v/c is the Mach number vector, and v is the fluid velocity vector 
x is the vector describing the ray path 

v is the unit vector tangent to the ray path, i. e. V = p/p and p = k/ko 
k is the wavenumber vector and k Q = uj/cq 

R is the vector describing the local geodesic elements of the wavefront, i. e. R 9 — (dx/dOo) ao and 
R a = (dx/da 0 ) Go 

Q are the conjugate elements of R 
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In Fig. 7 two typical examples of ray-tracing through a random field are shown. The upper graph 
presents an example of a random temperature field with a scale of L = 0.1 m. During the first 
correlation lengths the rays are only slightly distorted and the initial plane wave structure is still 
apparent. At greater propagation distances the deflection of the rays is more pronounced, 
neighboring rays begin to cross, strong concentrations occur at different lateral positions, and 
caustics appear. The lower graph shows the propagation through a random velocity field. The 
length scale is the same, and the index fluctuations for the scalar and vectorial fields are identical, 
i.e., e = —T'/2Tq = —v' x /co, where T' is the rms temperature fluctuation and v' x is the velocity 
fluctuation component in the propagation direction. The ray distortion for the velocity field is 
stronger than for the scalar field. The distance of caustics formation is shorter than for the scalar 
field. 


Ray tracing through a single realization of a 2D isotropic 


scalar field 


L = 0.1 m T'/2 T 0 = 5.882 10' 3 



x/L 


Single realization of a 2D isotropic vectorial field 


L = 0.1 m v' x / c 0 = 5.882 10' 3 



Figure 7: Examples of ray-tracing through a random temperature and velocity field. 
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NONLINEAR TRANSPORT EQUATION 


The main assumptions in the derivation of the nonlinear transport equation are: (1) only 
quasi-linear terms are retained in the derivation, (2) lossless propagation is assumed, and (3) the 
geometric acoustics approximation is applied. Using the assumptions, we derive a transport 
equation (Eq. 6), where s represents the arclength along the ray, f3 is the coefficient of nonlinearity, 
tl is the retarded time coordinate, and A is the cross-section of the ray tube. Through a 
transformation of the dependent and independent variable the transport equation takes the form of 
the Burgers equation (Eq. 9) for plane wave propagation in a homogeneous medium. The distortion 
distance variable Z describes the equivalent plane wave distortion for a wave propagating in a 
random medium. Equation 8 is integrated numerically along the rays and yields the equivalent 
distortion value. The Burgers equation (Eq. 9) is solved numerically by a Pestorius’ type algorithm 
[5]. The nonlinear distortion is applied in the time domain and an ad-hoc absorption is applied in 
the frequency domain. 


Nonlinear transport equation: 


Main assumptions: - quasi-linear 

- lossless 

- geometric acoustics approximation 


+ ( 6 ) 

Transformation of the dependent and independent variable 
Let Tl = Kp, where 


k = i*+mi m 

and 

f = (8) 

The transport equation is transformed into the homogeneous 
plane wave Burgers equation 


an _ go 

dZ p 0 cl 


nn t < = o 


0 ) 
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COMPUTATION STRATEGY 


The computation strategy involves the following steps. First, individual realizations of the turbulent 
field are contructed. Second, eigenrays that connect source and receiver are computed, and the 
properties K and Z are calculated along the eigenrays. Third, a Pestorius’ type algorithm is used to 
propagate a given waveform nonlinearly along the eigenrays. Finally, the waveforms associated with 
all eigenrays are combined to obtain the signal at the receiver. If a waveform passes through a 
caustic, the propagation proceeds up to the caustic, then a —90° phase shift is applied in the 
frequency domain, and propagation proceeds again. 


Strategy 


- create realizations of turbulent field 

- compute eigenrays and properties (K,Z) along the eigenrays 

- use Pestorius algorithm to propagate a waveform nonlinearly 
along the eigenrays 

- combine eigenrays to obtain signal at the receiver 

- if waveform passes through a caustic, propagation proceeds 
up to the caustic, then a —90° phase shift is applied 

in the frequency domain, and propagation continues 
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EXAMPLE OF EIGENRAY CALCULATION FOR ONE REALIZATION 


The upper graph in Fig. 8 represents the 7 eigenrays that connect the initial plane wavefront with 
the receiver located at 8 m on the x-axis. The second graph shows the values of K along the rays. 
When the eigenray passes through a caustic, the value of K drops to zero, i.e., an infinite pressure, 
a consequence of the linear geometric acoustics assumption. As is observed, one eigenray passed 
through two caustics, three rays passed through one caustic, and three rays did not yet pass through 
a caustic. The third graph presents the equivalent plane wave distortion distance Z. For the ray 
that passed through two caustics, the equivalent distortion distance is nearly 14 m when the actual 
propagation distance to the receiver is about 8 m. Thus, the nonlinear distortion along this ray is 
equivalent to that of a plane wave that propagated about 14 m in a homogeneous medium. For one 
ray that did not propagate through a caustic, the equivalent distortion distance is just slightly more 
than 2 m, and hence the nonlinear distortion will be very weak compared to the homogeneous case. 


Eigenrays 



Distortion distance variable 



Figure 8: Examples of eigenray and associated K and Z calculation. 
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PARAMETERS FOR THE NUMERICAL STUDY 


The parameters for the numerical experiment are shown in Fig. 9. Both temperature and velocity 
turbulent fields are constructed. Both fields have a length scale of L — 2.5 cm. The index of 
fluctuations e is identical for both fields. The rms velocity fluctuation is 2.5 m/s and the rms 
temperature fluctuation is 4.27 K. Each realization is constructed as a sum of 60 Fourier modes 
equally divided between a lower wave number value of 0.1 L and an upper value of 10 L. Statistics 
are calculated for an ensemble of 100 realizations. The incident acoustic wave is a plane N wave of 
peak pressure 500 Pa, duration 15 //s, and rise time 1 fis. The propagation medium is air, and 
classical thermo- viscous absorption is included in the Pestorius’ algorithm as well as 0 2 and N 2 
relaxation. 

Parameters for numerical study 


• Turbulence: 60 modes between 0.1 L and 10 L 

L = 2.5 cm 

v* = 2.5 m/s and T = 4.27 K 
100 realizations 

• plane N wave: 500 Pa peak pressure 

15 (is duration 
1 (is rise time 

• medium: classical absorption 

0 2 and N 2 relaxation 



Figure 9: Parameters for the numerical study. 
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AVERAGED RESULTS FOR THE PARAMETERS Z AND K 


In the upper graph of Fig. 10 the value of Z averaged over 100 realizations is shown. Three curves 
are presented, i.e., for the no turbulence, temperature turbulence, and velocity turbulence case. It is 
seen that the presence of turbulence always results in a lower value for the equivalent distortion 
distance. The effect is more pronounced when propagation distance increases, and stronger for the 
velocity turbulence. At a propagation distance of 2 m the equivalent distortion distance for the 
velocity turbulence is slightly more than 50% that of the no turbulence case. The lower graph 
presents the values of the amplitude factor K. The factor increases rapidly with propagation 
distance, and the effect is more pronounced for the velocity turbulence. 

Results for Z and K parameters 

Distortion distance averaged over 100 realizations 
2 

1.8 
1.6 
1.4 
1.2 

v (**0 | 

0.8 
0.6 
0.4 
0.2 
0 

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

distance source - receiver (m) 


Amplitude factor averaged over 100 realizations 
45 

40 

35 

30 

K 25 
20 

15 

10 

5 

0 

Figure 10: Averaged results for parameters Z and K . 
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EXAMPLE OF WAVEFORM CALCULATION 


An example of a waveform calculation is shown in Fig. 11. The total waveform at the receiver is 
shown in the upper graph. For this particular example 5 eigenrays exist. The lower graph presents 
the contribution of the waveform of each eigenray to the total signal. Two waves have passed 
through a caustic, and the other three have not passed through a caustic. The latter three arrive at 
the receiver at nearly the same time and are hard to distinguish. 


Example of calculation 



time (ps) 


Waveforms associated with each eigenray 



Figure 11: Example of waveform calculation. 
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WAVEFORM EXAMPLES 


In Fig. 12 several examples of waveforms calculated at the receiver at a distance of 18 or 20L are 
shown. The upper left graph presents the waveform at a distance of 18L in absence of turbulence. 
All other waves have propagated through a turbulent field. In general the distortion of the 
waveform is similar to that measured in spark-produced N waves (Fig. 3). 


Waveform examples 
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Figure 12: Waveform examples. 
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RESULTS FOR THE TOTAL WAVEFORM AT THE RECEIVER 


The average values of peak pressure and rise time of the total waveform calculated at different 
receiver distances are shown in Fig. 13. The rise time is calculated as the time portion between 10% 
and 90% of the peak pressure of the total waveform. According to this definition, we notice that the 
rise time is essentially determined by the differences in arrival time of the eigenrays when more than 
one eigenray is present (see, e.g., Fig. 11). The average is taken over 100 realizations. In the 
presence of turbulence, peak pressure is always decreased. There is no notable difference between 
the temperature and velocity turbulence. At a distance of 1 m a decrease of about 25% is observed. 
This result confirms the model experiment results [3] that peak pressure on average is always 
decreased by turbulence. The lower graph presents the average value of rise time. On average, rise 
time is always increased by turbulence. At a propagation distance of 1 m, a tenfold increase is 
observed. Again, this result confirms that of the model experiment. The effect is more pronounced 
for the velocity fields. The curves start to deviate from the no turbulence case when the rays pass 
through the first caustic. The shorter distance to the first caustic for the velocity fields explains the 
quicker departure from the no turbulence values. 

Results for total waveform at the receiver 


Average peak pressure for 100 realizations 



Average rise time for 100 realizations 

8 i * • • * 



0 i . . 1 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


distance source - receiver (m) 

Figure 13: Average values of peak pressure and rise time of the total waveform for 100 realizations. 
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HISTOGRAMS OF PEAK PRESSURE AND RISE TIME 


The histograms of peak pressure and rise time at different distances of propagation are shown in 
Fig. 14. The thick black vertical line in each graph represents the value for the plane wave case in 
absence of turbulence. The distribution for the peak pressure is asymmetric. Peak pressures smaller 
than the no turbulence value occur more often, but some large values of peak pressure are always 
present. The rise time histograms show that turbulence almost always causes the rise time to 
increase. When propagation distance increases, larger values for the rise time occur once the rays 
have passed through the first caustic (e.g., at a distance of 20 L). When rays pass through several 
caustics, rise times more than 20 times that of the no turbulence value are observed (e.g., at 40 L ). 


Histograms of peak pressure and rise time for 100 realizations 
of velocity turbulent fields 



Figure 14: Histograms of peak pressure and rise time for 100 realizations. 
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RESULTS FOR THE INDIVIDUAL WAVEFORMS ASSOCIATED WITH EACH EIGENRAY 


A second statistical calculation is performed. In order to filter out the effect of the difference in 
arrival time for all the eigenrays, we calculated the peak pressure and rise time of the individual 
waveforms associated with each eigenray. The waveform shown in Fig. 11 consists of 5 eigenrays. 
Instead of calculating the peak pressure and rise time of the global waveform, we calculate the peak 
pressure and rise time of the three waveforms that have not yet passed through a caustic. We do 
not take into account the rays that have passed through a caustic because the 10% to 90% criterion 
results in very large values. A better criterion for these waveforms would be a rise time value based 
on a maximum slope value. The result of the calculation is shown in Fig. 15. The average peak 
pressure is again always decreased by turbulence. Average peak pressure for velocity turbulence is 
lower than that for the temperature turbulence. The average rise time values are always increased by 
the turbulence. The difference is less marked when compared with the previous calculation. For the 
velocity turbulence a doubling of the rise time is observed at a propagation distance of 1 m. Again 
the effect is more pronounced for the velocity fields. The slope of the rise time versus distance curve 
is linear. The larger slope for the velocity fields is caused by the vectorial character of the fields. 

Results for individual contributions of each eigenray 


Average peak pressure for 100 realizations 




distance source - receiver (m) 

Figure 15: Average values of peak pressure and rise time of the waveforms associated with the 
eigenrays. 
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CONCLUSIONS 


We presented a model for the nonlinear propagation of sonic booms and spark-produced N waves 
through turbulence. A numerical experiment that simulates the propagation of spark-produced N 
waves through 2D temperature and velocity turbulent fields was performed. The results from the 
numerical experiment confirm the experimental observations of the model experiment in which 
spark-produced N waves propagated through a plane jet turbulence. An important conclusion is 
that the presence of turbulence on average reduces the nonlinear distortion. It is however possible 
that for a particular realization the nonlinear distortion is stronger. On average, peak pressure is 
always decreased by turbulence. Turbulence almost always causes the rise time to increase. The 
computed waveforms show similar distortion as that obtained in the model experiment. The effect 
of turbulent velocity fields is more pronounced than that of temperature fields and is caused by the 
difference in character, i.e., vectorial for the velocity fields versus scalar for the temperature fields. 
Future work includes a calculation with a more realistic turbulent energy spectrum and the 
incorporation of boundaries. 


Conclusions 


• Model presented for the nonlinear acoustic propagation 
through turbulence. 


• On average, turbulence reduces the nonlinear distortion. 


• On average, peak pressure decreases after propagation through 
turbulence. 


• Rise time is almost always increased by turbulence. 


• Effect of velocity turbulence is more pronounced than that of 
thermal turbulence. 

Figure 16: Conclusions. 
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A recently develped transonic small-disturbance model is used to analyze the interactions 
of random disturbances with a weak shock. The model equation has an extended form of 
the classic small-disturbance equation for unsteady transonic aerodynamics. It shows that 
diffraction effects, nonlinear steepening effects, focusing and caustic effects and random 
induced vorticity fluctuations interact simultaneously to determine the development of the 
shock wave in space and time and the pressure field behind it. A finite-difference 
algorithm to solve the mixed-type elliptic hyperbolic flows around the shock wave is 
presented. Numerical calculations of shock wave interactions with various deterministic 
vorticity and temperature disturbances result in complicate shock wave structures and 
describe peaked as well as rounded pressure signatures behind the shock front, as were 
recorded in experiments of sonic booms running through atmospheric turbulence. 

INTRODUCTION 

The review of various experimental and theoretical investigations of the interaction of 
shock waves and specifically sonic booms with free stream vortical or turbulent flows shows 
(Rusak and Cole 1 ) that this complex nonlinear interaction is still an open problem. 
Specifically, the improved simulation of sonic boom propagation through the real 
atmosphere requires a better understanding of the interaction of weak shocks with vortical 
perturbations and turbulence. 

Analysis of experimental data and theoretical approaches shows that in the case of 
the sonic boom, the shock waves near the ground are very weak, but still stronger than any 
acoustic wave. Also, flow fluctuations due to the atmospheric turbulence or vortical shear 
flows can become comparable to the shock weak strength such that locally the shock 
strength can be either strongly reduced or magnified and the shock wave front can be 
distorted significantly. Therefore, linearized acoustics and its second-order scattering 
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problem, or first-order linear theories of shock-vorticity interaction do not represent 
correctly the development of the weak shock and the pressure field behind it. However, in 
a coordinate system moving with the basic weak shock, the problem may fit the transonic 
framework. 

In a recent paper, Rusak and Cole 1 have presented a new extended transonic 
small-disturbance model to describe the interactions of random fluctuations with a weak 
shock wave. The model equation also has an extended form of the classic nonlinear 
acoustics equation that describes the propagation of sound beams with narrow angular 
spectrum (KKZ equation ) 2 " 3 and is similar to the model equation of Pierce 4 . The model 
shows that diffraction effects, nonlinear steepening effects, focusing and caustic effects and 
random induced vorticity fluctuations interact simultaneously to determine the 
development of the shock wave in space and time and the pressure field behind it. 

This paper summarizes the theory of Rusak and Cole 1 . A finite- difference algorithm 
to solve the mixed-type elliptic hyperbolic flows around the shock wave is also presented. 
Numerical calculations of weak shock wave interactions with deterministic vorticity and 
temperature disturbances describe both peaked or rounded pressure signatures as were 
recorded in experiments of sonic booms running through atmospheric turbulence 5 " 10 . 

A TRANSONIC SMALL DISTURBANCE MODEL 

The analysis of the linearized problem of the interaction of a weak shock with small 
disturbances shows 1 that it is an invalid approach when the flow perturbations are of the 
order of the shock strength. Therefore, a different approach has been developed to study 
the interaction of weak shocks with comparable random fluctuations in the flow (Rusak 
and Cole) L In a coordinate system moving with a basic given weak shock, the problem 
may fit the framework of transonic theory. A transonic small-disturbance model has been 
developed to analyze the flow across a basic weak shock running in the (— x) direction. 

A coordinate system attached to the basic shock is considered. The velocity vector (V), 

pressure (P), density ( p ) and vorticity (w) are described every where in the flow by: 
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V = U^, {i (1 + c 2/ 3 u + eui + e 4 / 3 U 2 + • • • ) 

+ j (evi 4 - e 4 / 3 v 2 + • • •) + k (twi + e 4 / 3 W2 + • • •)} 

( 1 ) 

P = p<x, (1 + e 2/3 P + epi + e 4/3 P 2 + • • •) 

P = Poo (1 + f 2/ 3 P + £p\ + C 4/ 3 p 2 + • • • ) 

U> = C 2 / 3 (ja/y + kw z ) + £ (iWxl + jWyl + kw z i) + • • • 


where = a 00 (l + K/2 e 2 / 3 ) is the speed of the basic shock (K > 0 ) and a^, p^, p w are 
the speed of sound, pressure and density of the unperturbed flow ahead of the shock. (e 2/ 3 ) 
represents the scale of strength of the basic weak shock where e < < 1. A rescaling of the 
x-coordinate and time (t) has also been considered: x* = x/e ^ 3 and t* = ta^e^ 3 , such 
that each of the terms in (1) is a function of (x*, y, z, t*). The rescaling in x means a 
stretching of the picture of the flow around the basic shock in order to capture the basic 
nonlinear effects that occur in the flow across the shock. The rescaling in time accounts for 
low-frequency unsteady perturbations in the flow. The constant K reflects that the speed 
of the basic shock wave is slightly higher than the speed of sound ahead of the shock. The 
substitution of Eqs. (1) into the continuity, momentum and energy equations results in 
(Rusak and Cole 1 ): 


u + P = f(y,z.t*) ( 2 a) 

^u + p = g(y,z,t*) (2b) 

2 lfiJ + (~ K - f +g -(7 + 1 ) u )^ + ^r + ^i l== -^^ 




(2c) 

£ 

II 

sia 5 

i 

dwi _ 1 dg 
dx * 7 dz 

(2d) 

i 

£ 

ii 

§> 

ii 

■ata> 

(2e) 
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where f and g are random induced fluctuations due to the free turbulence. The function g 
is related to the vorticity fluctuations in the flow and the function f is due to temperature 
or speed of sound fluctuations. Equations (2) show that the axial perturbation (u), 
pressure perturbation (p) and density perturbation (p), that are of the order of the shock 
strength (e 2 / 3 ), interact with the transverse velocity perturbations vi and wj, that are of a 
smaller scale (e). 

The substitution of u = g/7 + u in (2c), (2d) and (2e) results in a problem for solving 
a velocity potential function <p(x*, y, z, t*) where: 

u = 3^> y ' = 'By ’ W| = 3^’ p = -7 J$ r ( 3a ) 

2 <W + [K + * + f + (7 + 1) 4> xV - (* >y + * J — i (3b) 

In a conservative form Eq. (3b) is given by: 

[ 2 <P X * + ^ g] + + [( K + g/ 7 + f) <£ x * + (7 + 1) ^x*/ 2 ] x * ” (0y)y “ (0z)z = 0. (3c) 

The exact shock jump conditions (Ref. 11) must be satisfied along any shock surface 
x* - h(y, z, t*) = 0 that may appear in the solution. To the leading orders, they result in: 

[f] = 0, [g] = 0, (4a) 

~ 2 [*,*] ^ + ( K + 7 + f ) [K*\ + (7 + !) -y- + p y ] m ffi = °> ( 4b ) 

[M + Nf =0, [*,] + [**•] If =°- < 4c ) 

where [a] represents the jump across the shock property a, [a] = 7 b — a^. Equations (4a) 
show that to the leading order there is no jump in entropy across the shock, [S] = 0. 
Equations (1) and (2) also show that the local Mach number M/ at any point in the flow is 
given by: 


84 



( 5 ) 


M^-l = £ 2/3 u*, U* = j(7+ 1 )<^ X * + K +f + &}. 

The flow is locally supersonic when (7 4- 1) <f> * + K 4- f + g/7 > 0 , sonic when 
(7 + 1) 4- K + f + g/7 = 0, and subsonic when (7 + 1 ) ^ + K + f + g/7 < 0. 

Equations ( 3 ) and ( 4 ) are an extended version of the classic small-disturbance equation for 
unsteady transonic aerodynamics (Cole and Cook 12 ). The changes are due to the random 
terms g and f. Starting from given functions for f and g and initial conditions that describe 
a given basic shock, Eqs. ( 3 ) and ( 4 ) can be integrated in space and time to describe the 
development of the shock wave and pressure field behind it. 

An alternative approach has also been found by taking an x*-derivative of (2c) and 
using Eq. ( 2 a). The pressure perturbation (p) satisfies the equation: 


8 

dx* 



K ± f + g /7 _ 7 + 1 
2 2 7 p 




( 6 ) 


Equation (6) is an extended version of the classic KKZ equation that describes the 
propagation of nonlinear sound beams with narrow angular spectrum in an invisdd fluid 
(Zabolotskaya et al. 2 , Kuznetsov 3 ). Equation (6) also has a similar form to the model 
equation that has been recently developed by Pierce 4 using logical considerations only. 


Equations ( 3 ) and (6) show that diffraction effects, nonlinear steepening, focusing and 
caustic effects, and random induced fluctuations due to turbulence interact simultaneously 
to determine the development of the shock wave in space and time and the pressure field 
behind it. Turbulence tends to change the local speed of sound in the flow across the shock 
and through this effect to reduce or to magnify the strength of the jump along the basic 
shock (see Eq. ( 5 )) or to distort the shock front. These changes may result in unsteady 
motion of the shock front or in caustic vertices or in reflected shocks behind the incident 
wave that can produce the variety of pressure signatures of sonic booms that are measured 
in experiments 5 ' 10 . 
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FINITE DIFFERENCE SCHEME 


A finite difference algorithm to solve the unsteady mixed-type elliptic-hyperbolic 
flow around the shock wave has been developed. Murman and Cole 13 and Cole and Cook 12 
techniques are used. A fully conservative scheme that is based on the conservative form of 
Eq. (3c) is derived. In this way the difference equations also contain the shock relations 
(Eqs. (4)). 

Consider a uniform finite difference mesh (Ax*, Ay, Az, At*) in space and time, with 
points (x*, y, z, t*) labeled by (i, j, k, n). The results can be easily generalized to a 
variable mesh. Equation (3c) can be expressed in a conservative flux form for a box 
centered on a mesh point (i, j, k). Therefore, 


SF { ( 2 *>. + ¥) < ij.M, - ( 2 **' + 7 S ) , u, M -„ } 

+ A^{(( K + 7 g + f )^* + (7+1) ^* /2 )<‘* 

-((K + ig + f)^. + (7 + l)^/2) (i . i/4 j M 

_ ^z|(^ z ) — (0z) } ~ 0 • 

L ^ L l ( i,j,k + l/ 2,n) { i,j,k-l/ 2,n) J 


1/ 2,j,k,n) 


( 7 ) 


In Eq. (7), (<p y ) and (0 Z ) are always calculated from a centered expression. However, the 
approximation of (<p J strongly depends on whether locally, at a point, the flow is 
subsonic, supersonic, sonic or if it is a shock point. Extending References 12 and 13 
methodologies to our case and using Eq. (5), a centered approximation and a backward 
expression are given for u*: 
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u« =K + £(j,M) + ig(j,k,n) 

( i,j,k,n) 1 

+ Jjj* (*(i + l,j,k,n) - <p(i - 1 j,k,n)) 
u* b =K+f(j,k,n) + Ig(j,M) 

(i,j,k,n) ' 

+ ^pr(^( i ,j,k,n)-^( i - 2 ,j > k ,n)). (8) 

The local type of the flow is determined by the following table: 12 * 13 


condition 

u* c 

u* b 

local flow is 

1 

< 0 

< 0 

subsonic 

2 

> 0 

> 0 

supersonic 

3 

> 0 

< 0 

a sonic point 

4 

< 0 

> 0 

a shock point 


Table 1. Algorithm for local type of flow 


Equation (7) has been developed in a specific form according to the local type of the flow. 
The variety of difference forms for locally subsonic, supersonic, sonic or shock points are 
described in Rusak and Cole. 1 

Starting from initial conditions that describe a given shock wave in the space for 
t = 0 (or n = 0), and given temperature fluctuations f(y,z,t) and vorticity perturbations 
g(y,z,t), the various difference forms can be applied for n = 1 at any mesh point according 
to Table 1. They can be solved by an iterative point or line — or plane — relaxation 
algorithm until at any point max | G (i,j,k,l) | < 8 where 6 is a given small tolerance of 
convergence. Then 0 x *(i,j,k,l) can be calculated at any mesh point and the process is 
restarted for the next time step. In this way the shock motion and pressure field behind it 
can be integrated in space and time and the effect of various deterministic and random 
fluctuations f and g can be studied. 
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NUMERICAL RESULTS 


The finite difference algorithm to solve Eq. (3b) has been applied to a variety of 
two-dimensional and steady shock wave interactions with vorticity and temperature 
disturbances. Several problems have been studied where a nominal shock wave with 
K = 1.2 centered in the middle of the computational domain has been considered. The 
following boundary conditions were used: <f> + = 0 along inlet surface, <f>^ = — 1.0 along 
outlet surface (which satisfy the basic shock jump relations) and <j)^ = 0 along upper and 
lower surfaces. 


The first case considered vorticity fluctuations only where f = 0 and 


0 

g(y) = j 0.5 sin 47r(y - 1/4) 
0 


0 < y < 1/4 

1/4 < y < 3/4 

3/4 < y < 1 . 


( 9 ) 


Calculated pressure fields and profiles along various cross sections are shown in Figures 1 
and 2. The bending of the shocks is in phase with the velocity perturbations (i.e., a 
positive velocity perturbation produces a downstream deflection and visa versa) as was also 
described by Ribner 14 . Relative to nominal shock conditions (Figure 2a), the pressure 
jump decreases noticeably (approximately 40%) where the shock is bent upstream 
(Figure 2b), and increases significantly (approximately 40%) where the shock is distorted 
downstream (Figure 2d). The rippled wavefront leads to focusing and defocusing effects 
behind the shock, where diffraction effects also become dominant. 


In the second case only temperature disturbances have been considered where g = 0 

and 


r o 


%) = 


0.5 sin 167r(y — 1/4) 


l o 


0 < y < 1/4 
1/4 < y < 3/4 
3/4 < y < 1 . 


( 10 ) 


Pressure fields and profiles along cross sections are presented in Figures 3 and 4. Again, 
shock wave distortion results in significant local pressure jump reductions and increases 
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(see Figures 4b and c). Basic effects involved are similar to those described in the first 
case. Temperature or speed of sound disturbances strongly affect the basic shock as much 
as vorticity disturbances. 


The third case combined shock wave interaction with both vorticity and temperature 
disturbances as are given by Eqs. (9) and (10), respectively. Calculated results are shown 
in Figures 5 and 6. The shock front distortion is more pronounced and composed of the 
basic two harmonics of the imposed disturbances. It is found that shock pressure jump at 
certain locations is nearly eliminated due to the combined effects, resulting in a rounded 
pressure profile (Figure 6b). At other locations along the shock front, the pressure jump is 
significantly increased relative to the nominal shock jump with approximately 70 % (see 
Figure 6c). 


CONCLUSIONS 

A new transonic small-disturbance model has been developed where a rescaling of the 
axial coordinate and time has been considered to capture the basic nonlinear effects that 
occur in the flow across the shock. This model results in two alternative approaches: (1) an 
equation for solving a velocity potential function that is described by an extended version 
of the classic small-disturbance equation for unsteady transonic aerodynamics 12 , and (2) a 
nonlinear stochastic equation to describe the pressure field that is similar to the model 
equation recently presented by Pierce 4 using logical considerations only. This equation 
also has an extended form of the classic equation that describes the propagation of 
nonlinear sound beams with narrow angular spectrum 2 * 3 . 

Both approaches show that diffraction effects, nonlinear steepening, focusing and 
caustic effects and random induced turbulence fluctuations interact simultaneously to 
determine the development of a shock wave in space and time and the pressure field behind 
it. Turbulence fluctuations tend to change the local speed of sound in the flow across the 
shock and through this effect to reduce or magnify the strength of the basic shock. 

A finite difference scheme that uses Murman and Cole 13 finite-difference techniques 
for solving mixed-typed elliptic hyperbolic flows with shock waves has also been presented. 
Numerical simulations of two-dimensional and steady shock wave interactions with various 
deterministic vorticity and temperature disturbances have been shown. Results describe 



complicate shock wave structures and peaks as well as rounded local pressure signatures 
behind the distorted shock fronts. Similar signatures were recorded in the experiments of 
sonic booms running through atmospheric turbulence 5 ' 10 . 


ACKNOWLEDGEMENT 


This work was carried out with the support of the NASA Langley Research Center 
under Award NAG-1-1362 and NASA Training Grant NGT-S1113. The authors would 
like to acknowledge Dr. G.L. McAninch for funding and monitoring this research. 


REFERENCES 

1. Rusak, Z. and Cole, J.D., "Interaction of the Sonic Boom with Atmospheric 
Turbulence," NASA CP 10132, May 1993, pp. 65-91. (Also: AIAA Paper 
93-1943). 

2. Zabolotskaya, E.A. and Khorkhlov, R.V., Akust. Zhumal 15, 1, 1969, pp. 40. 

3. Kuznetsov, V.O., Akust. Zhumal 6, 4, 1970, pp. 548. 

4. Pierce, A.D., "Wave Equations and Computational Models for Sonic Boom 
Propagation Through a Turbulent Atmosphere," NASA CP-3172, 

February 1992, pp. 31-48. 

5. Hubbard, H.H., Maglieri, D.J., Huckel, V. and Hilton, D.A., "Ground 
Measurements of Sonic Boom Pressures for the Altitude Range of 10,000 to 
75,000 Feet," NASA TR R-198, July 1964. 

6. Maglieri, D.J., "Sonic Boom Flight Research — Some Effects of Airplane 
Operations and the Atmosphere on Sonic Boom Signatures," in Sonic Boom 
Research, Edited by Seebass, A.R., NASA SP-147, April 1967. 

7. Garrick, I.E. and Maglieri, D.J., "A Summary of Results on Sonic Boom 
Pressure- Signature Variations Associated with Atmospheric Conditions," 

NASA TN D-4588, May 1968. 

8. Downing, M.J., "Lateral Spread of Sonic Boom Measurements from U.S. Air 
Force Boomfile Flight Tests," NASA CP-3172, February 1992, pp 117-135. 

9. Willshire, Jr., W.L. and Devilbiss, D.W., "Preliminary Results from the White 
Sands Missile Range Sonic Boom Propagation Experiments," NASA CP-3172, 
February 1992, pp. 137-149. 

10. Lipkens, B. and Blackstock, T.T., "Model Experiment to Study the Effect of 
Turbulence on Risetime and Waveform of N Waves," NASA CP-3172, February 
1992, pp. 97-107. 


90 



11. Courant, R. and Friedrichs, K.O., Supersonic Flow and Shock Waves , Interscience 
Publishers, Inc., New York, pp. 297-302, 1948. 

12. Cole, J.D. and Cook, L.P., Transonic Aerodynamics, North-Holland, 1986. 

13. Murman, E.M. and Cole, J.D., "Calculation of Plane Steady Transonic Flows," 
AIAA Journal, Vol. 9, pp. 114-121. 

14. Ribner, H.S., "Cylindrical Sound Wave Generated by Shock-Vortex Interaction," 
AIAA Journal, Vol. 23, No. 11, 1985, pp. 1708-1715. 






perturbation pressure - p 
1.77682+ 

1 .59464 to 1.77682 
1.41245 to 1.59464 
1.23027 to 1.41245 
1.04809 to 1.23027 
0.865909 to 1 .04809 
0.683727 to 0.865909 
0.501545 to 0.683727 
0.319364 to 0.501545 
0.137182 to 0.319364 
ft -0.045 to 0.137182 


Figure 1. 


Pressure field for vorticity fluctuations (9). 
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Pressure field for combined vorticity and temperature disturbances. 
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INTRODUCTION 


Sonic booms received on the ground tend to be restricted to a region of finite lateral extent below 
the flight track. This occurs because of refraction and because the effective speed of sound, even with 
winds taken into account, decreases with altitude in the lower atmosphere. Not all rays proceeding 
initially downwards from the flight track within an allowable range of initial directions will reach the 
ground. The restricted region which can be reached by rays impacting the ground is known as the 
primary carpet. However, weak rumbles are heard in the nominal shadow zone beyond the edge of 
this carpet. A full wave theory is necessary for explaining waveforms in that region, and the present 
paper gives a matched asymptotic expansion technique for a suitable approximate full wave theory 
that involves a relatively small number of parameters. The outer solution is derived from the structure 
of the system of rays that impact near the corridor edge; the inner solution involves a solution of the 
parabolic equation and results in the special functions encountered in the diffraction of sound over the 
tops of hills. (Work supported by Armstrong Aerospace Medical Research Laboratory, Human Systems 
Division (AFSC), United States Air Force, Wright-Patterson AFB, Ohio 45433-6573.) 
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DETERMINATION OF ANGLE WITH CARPET EDGE 



i 



height above ground 


DETERMINATION OF EFFECTIVE SOUND SPEED 


r 


■\ 


Effective Sound Speed 



c(z) 

Ceff(z) = c(z) 



v(z) • egraze = horizontal wind component 
in direction of grazing ray 
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DETERMINATION OF RAY CURVATURE 
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THEORETICAL PREDICTION OF RAY TUBE AREA 
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EARTH-FLATTENING APPROXIMATION AND ITS CONVERSE 


Converse of Earth-flattening approximation: 


"Before" "After" 

c varies with z c independent of z 



• Use a set of curvilinear coordinates where grazing ray appears straight, 
but ground appears curved 

• Edge of primary carpet becomes line on top of cylinder 

• Grazing ray becomes horizontal straight line 

• Crosses the top of cylinder at oblique incidence with angle Xgraze 
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FOURIER TRANSFORM SOLUTION FOR WAVEFORM 
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FORMULATION BASED ON PLANE WAVE INCIDENCE 


Plane wave obliquely incident on a cylinder: 

• Trace velocity matching principle: 

P = pit- [ x/V f ],y,z ) 

a i a 

dx V f dt 

• Wave equation: 

d 2 p d 2 p /I 1 \ d 2 p 
dy 2 dz 2 \c 2 Vf) dt 2 

Results for normal incidence apply with transformation 

1 11 1. 2 
~ ? SU1 ^ 




WAVEFORM BEFORE ENTERING THE TRANSITION REGION 


~ N 

Waveform going into transition region: 


• Use N-wave of form 



P N = peak amplitude of N-wave 
T n = positive phase duration of N-wave 

• Determine Pn and T n by extrapolation of waveform received on the ground 
within the primary carpet, taking into account 

• Pressure doubling on reflection at ground 

• Geometrical spreading along ray tubes 

• Lengthening with propagation distance of due to non-linear effect 

• Extra decrease of Pn with propagation distance due to non-linear effect 
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FOURIER TRANSFORM SOLUTION FOR WAVEFORM 


Fourier transform solution: 


Let incident plane wave be of the form: 


Pine = Pn /n(^/Tn) 

T = t - [ y/c ] sin Xgraze - [x/c] COS Xgraze 

/ n(t/T n )= / e- lcor f N (coT N )T N dco 


Then the total solution is of the form 


r°° 

P = PnT n / z, a>)dco 

J —oo 


• where the transfer function TF satisfies 


9 2 TF 9 2 TF a) 2 . 2 _ 

T" H — H sm 2 y graze TF = 0 

9v 2 dz 2 c 2 g 
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TRANSFER FUNCTION FOR GROUND LEVEL SIGNATURES 


For points on the ground: 


• The transfer function is given by 


i[(o/c\y Sin Xgraze r( . n . 
TF(y, z,co) = e G(f, 0, q) 


/ OO 10 

.00 w[ (a) - 


qwi(a) 


wi(a) = 2n xll e inl ^ K\(ae ll7tiz ) 

H ~ y /[2 ^ <7 = ^cyl sin Xgraze) ^ pc/^s 


Ls.\ — 


Recall that 


»2/3 p2/3 

^cyl Xgraze Sin Xgraze 

(cue -1 sin Xgraze) 1/3 (CUC -1 ) 1 / 3 


2 

^cyl = Xgraze SU1 Xgraze 




WAVEFORM IN LIMIT OF HARD GROUND 


Hard ground limit: 


• Transmitted pressure reduces to 

/ OO 

r lW -i'' T « , A»)C«, o, o )dco 

- OO 


T—t- iy/c] sin Xgraze - [x/c] COS Xgraze 

* oo Ja% 


G($, 0,0) 


/ OO l 

4 

- 00^1 


(«) 


-da 


$ = ( j / i ™)«) i/3 

• Characteristic length for N-wave deterioration through and beyond the tran- 
sition region 

f^tran = S in Xgraze ^ 













To study deterioration of front of shock, 

sufficient to take incident wave as a step function: 


P 1 

p = — — Im[e -i<WT G(£, 0, 0)]d 

n Jo o) L 


o>--G(0,0, 0 )P N 


• To examine the contribution of the higher-frequency components, it suffices 
to use the creeping wave approximation: 

„ , e~ i7l/6 K% 


G 1.8325e " j 

where K = 1.0188 . Result is (for y > 0) 


as £ — > +00 


P — P N ^ creep (T / ^char) ^ ^ (®* 

<tW») = — r ^- [V3/2ln ' ,! sin(a„ - [l/ZlO^yo 

71 Jq 

(1.088) 3 y 3 

^char — ~ tt; 

2 cR* sm X graze 





HINTS FOR ANALYSIS OF FIELD DATA 


For analysis of field data: 

• Need at outset to identify 

• nominal location of carpet edge, origin of y 

• Mach number relative to sound speed at ground, to get x% raze 

• Grazing ray radius of curvature R graze, from meteorological profiles 

• N-wave amplitude Pm before reflection at ground 

• Positive phase duration Tm at the carpet edge 

• Theory predicts that peak waveform versus extra lateral distance should fall 
on a “universal curve” if one plots p^/Pn versus y/Ltran 

• Theory predicts that waveform shape is similar for all waveforms recorded 
at same value of y/L u ^ 




HOW WAVEFORMS EVOLVE BEYOND CARPET EDGE 


Rough characterization of wave evolution beyond the transition region. 


1 


Begin with an N-wave of amplitude of 1.399 Pn at the carpet edge (y = 0). 


|G(£, 0, 0)| drops to 1/2-th of its zero-frequency value at f = 1.05. 


Frequencies which propagate with negligible attenuation to lateral distance 
y from carpet edge are those for which 


co < 


^ C ^graze si** 3 Xgraze 


• Propagation to lateral distance y from carpet edge is roughly equivalent to 
passing an N-wave of peak overpressure 1399 Pn through a low-pass filter 
with cut-off frequency 

2c/?^raze S ^ 3 Xgraze 
^cut-off = 5 
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Accurate prediction of human response to sonic booms from proposed HSCT aircraft depends 
on a knowledge of the waveshape and risetime of the boom at the ground. In previous work, we have 
developed a numerical technique to predict the combined effects of molecular absorption and finite 
wave distortion on the sonic boom as it propagates from the aircraft to the top of the turbulent boundary 
layer. We have more recently developed a scattering center based model to calculate the effects of tur- 
bulence on the sonic boom waveform as it propagates through this boundary layer. Calculations have 
been performed using single scales of turbulence and compared to measurements at Edwards AFB in 
the late 1960s. A model of the atmosphere involving two scales each for convective and mechanical 
turbulence has been developed and fit to meteorological data collected during JAPE 2. Scattering calcu- 
lations employing this model underpredict the number of unperturbed waveforms. In order to develop 
a more realistic model of the atmosphere, the JAPE 2 meteorological data has been fit to a von Karman 
spectrum. Results of scattering using this multi-scale model will be presented. The combination of 
finite wave effects with turbulent scattering predictions includes the principal effects of the atmosphere 
on the sonic boom from the HSCT. 
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INTRODUCTION 


The prediction of the average environmental impact of the HSCT requires accurate modeling of 
the processes affecting the sonic boom waveform and risetime. We have used the enhanced Anderson 
algorithm to predict the risetime and waveshape of sonic booms under non-turbulent conditions. This 
method can also be used to predict the risetime and waveshape at the top of the turbulent planetary 
boundary layer. 

The enhanced Anderson algorithm includes all finite wave effects and the vibrational relaxation 
effects of N2, 02, and CO2 in combination with atmospheric H2O. This algorithm has been compared 
to data from explosions 1 and sonic booms 2 and has been tested against measurements of high intensity 
ballistic waves from rifles and from tank guns 3 . In addition, the results of this calculation for quasi 
steady shocks agree with the results from the enhanced Burgers’ Equation 4 ’ 5 . 

Figure 1 presents the results of the application of the enhanced Anderson algorithm to a 
predicted HSCT waveform 6 . We emphasize that the key parameter in determining the risetime of the 
sonic boom is the absolute humidity. 

Under turbulent conditions, the risetimes of sonic booms are scattered and are occasionally as 
large as ten times the risetimes calculated from vibrational relaxation considerations. It is clear that 
turbulence is the cause of the increased risetime and peculiar waveforms observed. Analytic techniques 
have been used to estimate the increase in average risetimes 7 ’ 8 - 9 and to calculate perturbed waveforms 
due to focusing and defocusing of the waves by turbulence 10 . In such calculations, it is usually 
necessary to assume a single strength and turbulence scale representative of the atmospheric turbulence. 
The largest turbulence effects are usually identified when the largest scales are chosen as typical. 

We have chosen a different approach to calculating the effects of turbulence on sonic boom 
risetimes and waveforms based on a simple scattering center-based theory. The scattering center-based 
method accurately predicted the effects of turbulence on the coherence of continuous wave signals 
above natural ground surfaces 11 . 


METHOD 

The scattering center-based technique resolves atmospheric fluctuations into a sum of discrete 
spherically symmetric Gaussian “turbules”. The total effect of the atmosphere is then calculated by 
summing up the scattering amplitudes. See Figure 2 . The scattered amplitudes are calculated using the 
first Born approximation. If the complex pressure at the receiver is written as: 
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( 1 ) 


N 

P B (?) = p 0 (?) + ^ \ff 

i = i 

where the superscript B refers to the first Bom approximation, p 0 (?) is the unperturbed spherical wave, 
and N is the number of turbules, then 


e ik(r M +r tt ) 

l 

r st + r tr 

1 - ia. 


»-Ck?s 2 /4 


( 2 ) 


where 


C = (1 - cos 0 O ) 2 + sin 2 0 O 


1 

1 - ia 


(3) 


and 



_L 

r s t 



(4) 


s defines the 1/e 2 contour of the turbule, qi is the index of refraction profile strength, and 0 O is the 
scattering angle. The geometry is indicated in Figure 3. 

The initial research on continuous wave propagation modeled the atmosphere as a random sum 
of identical turbules. This single scale calculation was extended to impulse propagation with promising 
results. 12 The impulse is Fourier transformed into the frequency domain and the total scattered compo- 
nent at each frequency is calculated. Then the inverse Fourier transform yields the time domain wave- 
form. The single scale calculation (s = 10m, 30m or 100m) with a fluctuating index of refraction of 
<|i 2 > = 10 x 10' 6 predicted spiked and rounded waveforms and predicted risetimes as large as 10 ms. 
These results encouraged us to analyze the results of the JAPE-2 tests 13 ’ 14 using the scattering center- 
based model. 


ANALYSIS OF JAPE-2 DATA 

The JAPE-2 tests consisted of simultaneous measurement of sonic boom characteristics and 
meteorological measurements. The wind and temperature fluctuations were measured at heights up to 
30m using sonic anemometers and hot wire anemometers. The sonic boom data was analyzed by 
Wiltshire, Garber and DeVilbiss 14 and provided as computer files. The turbulence data was analyzed 
by Bass, Boulanger, Olsen and Chintawongvanich 15 . 
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a.) Two Scale Model 


Examination of the data showed that a single scale model of the atmosphere could not fully 
describe the turbulence above the ground. The time correlation of the fluctuation quantities was fit to a 
two scale model. See Figure 4. Table I displays the results of the analysis for a moderately turbulent 
day during JAPE-2. 


Table I. Example of the Two Scale Model Applied to Atmospheric Data 



Wind driven 1 

Wind driven 2 

Temp, driven 1 

Temp, driven 2 

Size (meters) 

117 

11 

74 

8 

Number of 
Eddies 

1 ; 

90 

1 

233 

V 2 ) 

0.54 10- 5 

0.25 10- 5 

0.5 10- 6 

0.4 10- 6 


The scattering calculation was performed by summing the results of four calculations - one for 
each scale size. The input waveform to the scattering calculation was an N-wave propagated from the 
flight altitude to the top of the turbulent layer using the enhanced Anderson algorithm. The results of 
the Anderson algorithm agree moderately with the measurements taken under low turbulence condi- 
tions. See Table II. 

Table II. Comparison of Measured and Predicted Waveform Parameters for the T-38 



Measurements for the low 
turbulence case 

Calculations using the 
Anderson algorithm 

Peak overpressure (psf) 

0.71 

0.88 

Risetime (ms) 

0.32 

0.33 


Figure 5 compares the results of the measurement and prediction for T-38 overflights under 
moderate turbulence conditions. Although the scattering center model produces a wide distribution of 
risetimes, it does not predict the shift of the histogram maximum to 2 ms; rather the maximum remains 
at the unperturbed value of 0.3 ms. It is believed that this is due to the use of two relatively large scales 
to represent the atmospheric turbulence. The scattering from large turbules is predominantly in the 
forward direction, and large turbules are relatively sparse, so that it is easy to “miss” the receiver with 
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the scattered wave. The four scale model does, however, represent a significant improvement over the 
single scale model. 

b.) von Karman Spectrum Model 

The fit of the autocorrelation to two scales rather than one improved the prediction of risetimes 
significantly. The high occurrence of unperturbed risetimes indicated that smaller and intermediate 
scales were needed to fully describe the scattering of sonic booms by turbulence. 

De Wolf 16 presented a technique for simulation of a turbulent atmosphere obeying the von 
Karman spectrum in terms of the number density of turbules. 

The general form of a 3-D von Karman spectrum is given in terms of frequency by: 

♦ (f)- a£ -rr (5) 

where: 

/ \ 2/3 

a = 4^C?(^) (6) 


and 



(7) 


The coefficients a and b are determined by fitting a function (j>(f) through the measured spectra. See 
Figure 6. 

The fit parameters are then used to determine n (s), the number density of turbules of size s 
needed to model the fluctuating atmosphere. De Wolf’s model was originally developed to predict 
second moments of a scattered field and therefore is designed to reproduce only second moments of the 
fluctuation fields. Higher moments must be accurately represented to express the temporal characteris- 
tics of an impulse. De Wolf used an index of refraction maximum for each turbule of ± 1.0 and em- 
ployed a very sparse distribution. We have varied the product of qi 2 and n (s) until the model distribu- 
tion approximates the measured second and fourth moments <|J, 2 > = 9.6 x 10' 6 , <|X 4 > = 2.5 x 10 10 . 
The variation of calculated <|i 2 > with number of turbules and qi 2 is shown in Table HI. 
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Table DL Calculate <ji 2 > and <|i 4 > as a Function of Number of Turbules 


Number of 
Turbules 

Percentage of 
Volume 

qi 2 

■EM 

V 4 ) 

42000 

8 

1.5 10- 4 

1.2 10- 5 

1.0 10-9 

63000 

12 

1.0 io- 4 

HBEKSflli 


95000 

18 

6.7 10- 5 

1.0 10-5 

6.9 io- 10 

127000 

24 

5.0 10- 5 

1.1 10-5 

5.5 10- 10 

254000 

48 

2.5 10- 5 

9.9 10- 6 " 

■■EE1ES9HH 


The turbule spatial and size distribution for each realization is determined by Monte Carlo 
methods. The index of refraction fluctuations along a straight line has been compared to the corre- 
sponding measured values and exhibits similar fluctuation scales and displacement 

The second improvement to the scheme was the use of the measured height of the Planetary 
Boundary Layer in the calculation. Figure 7 displays the temperature versus height curve for one flight 
during JAPE-2. One sounding is taken with the tethesonde going up and the inversion height is 400m, 
the other trace is the tethesonde coming down 30 minutes later and the inversion height is at 670m. The 
turbulent layer thickness at the time of the later sonic boom measurement was extrapolated from this as 
750m. 

The results of this calculation for 20 realizations are displayed in Figure 8. The maximum 
occurrence risetime shows a shift away from the non-turbulent risetime of 0.3 ms. The smaller and 
intermediate scales of turbulence have a significant effect on the risetimes of sonic booms. It is clear, 
however, that the shift is not large enough to match the measured data in Figure 5a. 


CONCLUSION 


The enhanced Anderson algorithm provides a good prediction of waveshape and risetime of the 
HSCT at the top of the Planetary Boundary Layer. 

The scattering center-based model can be extended to predict distorted wave shapes and longer 
risetimes. At this stage, the scattering based model does not predict long enough average risetimes, but 
does show that smaller and intermediate scales are important in increasing the average risetimes. 

The larger scales are the source of the dramatically distorted waveforms, but are not the source 
of the shift in average risetimes. The scattering center-based calculation allows the quantitative investi- 
gation and modeling of the turbulence effects discussed qualitatively by Crow, Plotkin and George, and 
Pierce. 
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Pressure (atm/500) 



Figure 1. a.) Predicted HSCT waveform close to the aircraft. 
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Autocorrelation 



33 







Percentage of Occurrence Percentage of Occurrence 


30 


20 



Rise Time (ms) 

Figure 5. a.) Measured risetime distribution for the T-38 aircraft. 



b.) Predicted distribution using the two scale model. 
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As U.S. aircraft manufacturers now have focused their HSCT efforts on overwater super- 
sonic flight, a great deal more must be known about sonic booms propagating overwater and in- 
teracting with the ocean. For example, it is thought that atmospheric turbulence effects are often 
much less severe over water than over land. Another important aspect of the overwater flight 
problem is the penetration of the sonic boom noise into the ocean, where there could be an envi- 
ronmental impact on sea life. 

T his talk will present a brief review on the penetration of sonic boom noise into a large 
body of water with a flat surface. It has been determined recently that faster supersonic speeds 
imply greater penetration of sonic boom noise into the ocean. The new theory is derived from 
the original Sawyers paper and from the knowledge that for level flight a boom’s duration is pro- 
portional to the quantity M/(M 2 — l) 3 / 8 where M is the Mach number. It is found that for depths 
of 10 m or less, the peak SPL varies less than 6 dB over a wide range of M. For greater depths, 
100 m for example, increased Mach numbers may increase the SPL by 15 dB or more. 
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INTRODUCTION 


Times are certainly changing in atmospheric propagation research for the NASA High 
Speed Research program. During the last few years research has focused on the effects of at- 
mospheric turbulence on the propagation and distortion of sonic boom waves. This research has 
yielded many important results. 

However, the U.S. aircraft manufacturers who are preparing for the construction of super- 
sonic commercial aircraft have shifted their emphasis away from overland flight. They are now 
centering on flying overwater to essentially eliminate the sonic boom impact on people. 

This decision has important repercussions for propagation research. It is known that atmo- 
spheric turbulence effects are much less severe over water than over land. Over land the ground 
heats up during the day as the sun shines, creating thermals which greatly increase atmospheric 
turbulence. The ocean surface, in comparison, does not heat up. The ocean acts as a great heat 
reservoir, and the water churning over and over keeps the surface temperature nearly constant. 
No thermals affecting turbulence are created. 


Changing Times 


• Last few years, HSR focusing on atmospheric turbulence. 

• U.S. aircraft manufacturers say 

- no to overland flight 

- yes to overwater flight. 

• Atmospheric turbulence much less severe over water than over land. 





NEW TECHNICAL CHALLENGES 


Since atmospheric turbulence does not seem to be as pressing an issue for the overwater 
program as it did for the overland program, one may ask: What are the technical challenges that 
will be encountered in overwater flight? 

Just as in the planning of any large project, one must be concerned with environmental im- 
pact. Thus the environmental noise impact of sonic booms on wildlife should be addressed. It 
turns out that the local wildlife for overwater flight consists of sea life, and particularly marine 
mammals such as whales. 

All whales breathe air, and must surface periodically. Most species spend the great majority 
of their time within 100 m or so of the ocean surface, many within the top 25 m. 

One may then ask: How much sound gets from the air into the water? What would one 
hear if just under the water’s surface? 


Overwater brings new challenges 


• Look at environmental noise impact on wildlife. 
Overwater local wildlife: marine mammals 



• How much sound gets from the air into the water? 





SOUND PENETRATION THEORY 


It turns out that one can make predictions with reasonable certainty, as the acoustical theory 
has been well understood for many years. One first notices that the characteristic impedances for 
water and air differ greatly, approximately 415 for air and 1,500,000 for water. This means that 
a plane sound wave propagating in air directly toward the surface of the water at an angle of 0 
degrees, called normal incidence, would have 99.8% of its energy reflected. Very little propagat- 
ing sound energy gets into water. 

Further it turns out that the sound from a sonic boom is not normally incident but is inci- 
dent at an angle greater than the critical angle, which for the air - water interface is 13.2 de- 
grees. This means that 100% of the incident energy is reflected. One might think that this means 
no sound gets into the water, but this is not the case. 

When a force is applied normally to the surface of some material, the material must push 
back or the surface will be moved. So obviously when a sonic boom is incident on the water 


Brief Review of Sound Penetration Theory 

• Characteristic impedances differ greatly: 

air: c\p\ = 415 kg/(s m 2 ) 
water: C 2 P 2 = 1,500,000 kg/(s m 2 ) 

• For typical HSCT flight, sound will be incident at angles greater than 
the critical angle, 13.2°. 

==4> all sound power reflected 

=4* but must still match pressure boundary condition 
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surface, the water must push back on the surface. What this means is that the acoustic pressure 
on both sides on the ocean surface must be the same. This boundary condition is important be- 
cause it implies that there will be substantial sound in the water although all of the energy is re- 
flected. 

If one uses typical numbers for a projected high speed civil transport, such as Mach num- 
ber of 2.1 and boom duration of 0.3 second, one can straightforwardly predict what the acoustic 
pressures will be under the surface of the water. Mach 2.1 implies that the incident angle will be 
28.4 degrees from nonpal, which implies that there will be a pressure wave in the water whose 
amplitude decays with depth. The “fundamental” frequency of a 0.3 second boom is 3.33 Hz, 
and the decay of this frequency is indicated in the figure. As one can see, given a relative am- 
plitude of 1 at the ocean surface, the amplitude decays slowly with depth. There is significant 
noise penetration beyond a depth of 50 m. Of course a sonic boom is composed of many com- 
ponent frequencies. Higher frequencies will penetrate less far beneath the surface, while lower 
frequencies will penetrate further. 


Matching the boundary condition with typical HSCT numbers: 

• Speed: Mach number, M = 2.1 

angle of incidence = sin -1 (1/M) = 28.44° 

=> exponentially decaying pressure in water, e ~ A2n f z 


• Duration, T = 0.3 s 

=4> “fundamental” frequency, /o = 3.33 Hz: 
Much faster decay for higher /. 


relative amplitude 


depth, m 
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By performing a superposition of component frequencies of a real sonic boom, one can fur- 
ther predict what the received noise would look like at various depths under the ocean surface. 
For example, using the 0.3 second duration boom incident at Mach 2.1, the Fourier analysis re- 
sults are shown below in the figure. It is assumed that there is a perfectly N shaped sonic boom 
incident on the ocean surface. One can see that at 2.5 meters depth, the tips of the N are slightly 
rounded, but the waveform is mostly preserved. At 25 meters depth the waveform’s peak ampli- 
tude is decreased to about 45% of its amplitude at the surface, and the waveform is rounded. At 
a depth of 120 meters or more, the waveform is greatly smoothed but still has an amplitude of as 
much as 15% of its peak at the surface. 


But a sonic boom is a pulse (Mach 2.1): 



time 
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In 1968 the first theory to account for this penetration of sonic boom noise was developed 
by Sawyers. His theory is in terms of several nondimensionalized variables to represent time, 
distance, and depth. The theory assumes that the ocean is flat and is so deep that one need not 
account for bottom reflections. The assumption is also made that N wave shaped booms are in- 
cident on the ocean surface. 


Sawyers’ Theory 


Boom penetration theory due to Sawyers (1968) : 


i fr + £- 1 


^surface 


= (2 r + 2£ — 1) tan' 


/o i ot 1W -i { r + , /-l r c 2 + ( r + £) 2 

~(2r + ae-l) tan ^ + (t + J _ ip 

where r, and ( are nondimensionalized t, x, and z. 


• Theory assumes ocean is perfectly flat, ocean is very deep, and 
booms are perfectly N shaped. 
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In 1970 Cook also closely examined the theoretical aspects of sonic boom penetration into 
the sea. Cook suggested several minor improvements to Sawyers’ theory, but Cook noticed that 
the small differences between the theories were unlikely to be seen in field experiments. 

In the early 1970s two separate laboratory experiments validated Sawyers’ theory. Waters 
and Glass exploded small charges in air, modeling sonic booms, above a pond of water and mea- 
sured the response of hydrophones below the water surface. Their measured waveforms closely 
matched those predicted by Sawyers’ theory, even though their incident waves were small explo- 
sions instead of sonic booms. 

Further Intrieri and Malcolm investigated the waveforms created by small supersonic pro- 
jectiles. They sent such projectiles through the air over an aquarium and measured the acoustic 
pressure waveforms in the water. Intrieri and Malcolm’s results were also in good agreement 
with Sawyers’ theory. 


• Sawyers’ theory was later examined by Cook (1970), who suggested 
minor improvements. 

• Theory validated by laboratory experiments of Waters and Glass 
(1970, 1972) and Intrieri and Malcolm (1973). 
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TODAY’S TALK: EFFECT OF AIRCRAFT SPEED 


What hasn’t been completely understood from Sawyers’ theory alone is the effect of aircraft 
speed on the penetration of sonic boom noise into the ocean. As one examines Sawyers’ nondi- 
mensionalized variables one can see that V, the aircraft speed, and T, the duration of the sonic 
boom, appear several times. Clearly to apply the Sawyers’ theory correctly, one must have accu- 
rate speed and duration information. 


Effect of aircraft speed? 

(The major topic of this talk.) 

• Examine Sawyers' nondimensionalized variables: 

C = z/(mT) r = t/T £ = x/(TV) 

where 

m = V (1 — V 2 /^ ater ) _1/2 , 

V - aircraft speed, and T - boom duration. 

• Call p = p I ^surface ■ 
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It turns out that the duration of the boom, T, is a function of the aircraft speed, V. This re- 
lationship is also more complicated than one might expect. Linear theory would predict that 
T=L/V where L is the length of the aircraft. Hence for a fixed velocity, longer aircraft (spatially) 
have longer sonic booms (temporally). 

However, such linear theory neglects nonlinear acoustics effects. Because the sound from 
a sonic boom is reasonably loud, particularly near the aircraft, the sound waveform of a sonic 
boom lengthens as it propagates toward the ground. Thus one must account for this finite ampli- 
tude acoustical nonlinearity when determining the correct relationship between T and V. 





Several authors have shown that the duration of a sonic boom due to a projectile has the 
form shown below, as a function of Mach number, which is equal to V divided by the speed 
of sound. This T versus V relationship assumes that the aircraft is flying at a fixed altitude in 
steady flight. The theory has been compared to experimental data with good success. In spite 
of this, the theory assumes that there are no lift effects as one has for an airplane. More compli- 
cated theories are available. 


Finite amplitude nonlinear acoustics theory predicts: 

M 

T ~ K {M 2 - 1)3/8 

where k is a constant; Pierce (1989), Maglieri and Plotkin (1991). 
Assumptions: 

- level flight at fixed altitude 

- lift effects ignored 

This theory is pretty good in describing T = f(V). 
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Given that we now have a model relationship between T and V, and know how these factors 
influence the penetration of sonic boom noises into the ocean, it is time to determine what will 
be the maximum acoustic pressures which occur under the water’s surface. One can determine 
the maximum pressure by taking the time derivative of the acoustic pressure from Sawyers’ the- 
ory using the scaled time variable, setting this result to zero, and determining which value of 
scaled time satisfies the resulting equation. This procedure is necessary because the value of 
scaled time at which the maximum pressure occurs varies with depth. 



Call this scaled time tmax(C)- 
The maximum value of p is then 


Pm ax = p(£ = 0, tmax(C)) 


Shown in the figure below is the maximum pressure plotted as a function of Mach number 
for three different depths. The depths are 1 m, 10 m, and 100 m below the ocean surface. These 
curves are based on some SR-7 1 aircraft sonic boom experimental data in the possession of the 
author. The curves are .based on an SR-71 traveling at a Mach number of 2.6 producing a nearly 
perfectly shaped N wave having a duration of 200 milliseconds. 

As one can see the maximum scaled pressures generally increase with increased Mach num- 
bers. The effect is more noticeable at deeper depths. Along the surface of the ocean the scaled 
maximum acoustic pressure would be 1 for all Mach numbers. 




Shown in this figure is the same information, but on a relative decibel scale. The curve for 
each depth is scaled by the maximum pressure at Mach 2.6. One can see that for depths of 10 m 
or less, the maximum sound pressure level, SPL, will vary less than 6 dB over a wide range of 
Mach numbers. For greater depths, 100 m for example, increased Mach numbers may increase 
the SPL by 15 dB or more. 



Instead of peak levels, one may also be interested in the time domain representations of the 
boom noise. Shown below are waveforms of scaled acoustic pressure versus scaled acoustic 
time. At the surface of the ocean the scaled pressure would have a value of unity and the dura- 
tion of the boom in scaled time units would also be unity. The three columns represent an air- 
craft having speeds of Mach 1.4, 2.4, and 3.4. In each column the waveform is shown at depths 
of 1 m, 10 m, and 100 m. 

One can clearly see that the waveforms at a depth of 1 m are nearly N waves. Further at 
10 m the waveforms are similarly shaped and have amplitudes well in excess of 50% of the cor- 
responding amplitude at the surface. The waveforms at 100 m are significantly smoothed out, 
except for high Mach numbers where a sizable waveform still exists. 





CONCLUSIONS 


The major conclusion from this study is that faster flying supersonic aircraft produce sonic 
booms which penetrate more deeply into the ocean. The acoustic pressures experienced under 
the ocean’s surface increase with increased Mach number. 

The key to making these predictions using Sawyers’ theory is an accurate functional rela- 
tionship between the boom duration, T, and the aircraft velocity, V. Although the functional re- 
lationship used here has wide agreement with experimental data, it certainly could be improved 
upon to account for other factors such as aircraft lift and geometry. 

This research is just the beginning, however. Many other factors must be taken into account 
to ascertain the sonic boom noise impact/non-impact on marine mammals. First of all, Sawyers’ 
theory assumes that perfectly N shaped waves are incident. A similar theory should be devel- 
oped for more realistically shaped sonic boom waveforms. 


Conclusions 

• Faster flying aircraft ==> increased penetration of boom noise. 

• Using the correct T = f(V ) relationship is the key. 

• More work needs to be done; this is only the beginning: 

Issues: Non N shaped waves. Include lift effects. 





AN IMPORTANT UNADDRESSED ISSUE 


In addition, another important issue should be addressed. The real ocean surface is rarely, if 
ever, flat. One should develop a boom penetration theory which takes into account the curvature 
of the ocean’s surface. 

This curvature causes the rippling surface of the ocean to act as a series of lenses, focusing 
and defocusing the sound under the surface. The focus spots can be considered acoustically hot, 
where the peak sound pressure levels could be substantially higher than in other regions below 
the surface. There has been no published research regarding this focusing and defocusing effect. 


Another issue: 

How often is the ocean perfectly flat? 

Instead have focusing and defocusing for real sea surface: 


incident boom, 
propagating 



\i / i \\f ! 

COLD 'HOT 'coLD HOT' 



It would be standard to represent the sea surface as a superposition of Fourier modes as a 
first pass model. Then one could investigate the focusing effect by a number of analytical and 
numerical techniques. 

For weak waves on the sea, called wind waves, analytical solutions would be tractable via 
perturbation expansion techniques available in modern symbolic algebra packages. 

For wind waves of higher amplitude one could look into either analytical or analyti- 
cal/numerical solutions to the Kirchhoff-Helmholtz integral equations. Such integral equations 
have been used to examine the propagation of constant frequency sound from air into water with 
a rippled surface. The extension of this methodology to sonic booms seems straightforward. 

Further one can also perform finite difference calculations with conformal grids to model 
the penetration of sound in the air into the ocean. Such calculations would be invaluable for 
validating the perturbation expansion and Kirchhoff approaches. Hence, there are a number of 
methods available today for predicting the focusing effect. 


Describe sea surface by Fourier superposition of many modes. 

How to handle: 

• Weak wind waves: Analytical perturbation solution via Mathematica. 

• Higher wind waves: Kirchhoff-Helmholtz integral formulation. 

(Others have looked at this but not for sonic booms.) 

• For comparison: Finite difference time domain calculations. 
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INTRODUCTION 


Sonic-boom propagation from flight level to ground is influenced by wind and speed-of-sound 
variations resulting from temperature changes in both the mean atmospheric structure and small-scale 
perturbations (refs. 1-5). Meteorological behavior generally produces complex combinations of 
atmospheric perturbations in the form of turbulence, wind shears, up- and downdrafts, and various wave 
behaviors. Differences between the speed of sound at the ground and at flight level will influence the 
threshold flight Mach number for which the sonic boom first reaches the ground as well as the width of the 
resulting sonic-boom carpet (refs. 6-8). Mean atmospheric temperature and wind structure as a function of 
altitude vary with location and time of year. These average properties of the atmosphere are well- 
documented (for example, refs. 9 and 10) and have been used in many sonic-boom propagation 
assessments. In contrast, smaller scale atmospheric perturbations are also known to modulate the shape 
and amplitude of sonic-boom signatures reaching the ground, but specific perturbation models have not 
been established for evaluating their effects on sonic-boom propagation. The purpose of this paper is to 
present simple examples of atmospheric vertical temperature gradients, wind shears, and wave motions 
that can guide preliminary assessments of nonturbulent atmospheric perturbation effects on sonic-boom 
propagation to the ground. The use of simple discrete atmospheric perturbation structures can facilitate the 
interpretation of the resulting sonic-boom propagation anomalies as well as intercomparisons among varied 
flight conditions and propagation models. 


Outline 

• Introduction 

• Example atmospheric profile 
• Wind shear probabilities 

• Temperature gradient statistics 
• Single mode wave model 
• Closing remarks 
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INTRODUCTION (cont'd.) 


An observed atmospheric profile is used to illustrate discrete atmospheric layers with strong wind 
shear and high-temperature gradient values. Statistics for vertical wind shears measured by FPS-16 
radar tracking Jimsphere balloons and by the conventional upper-air rawinsonde measurements (refs. 
10-14) are presented. These wind statistics and data for vertical temperature gradients in the lower 
stratosphere also indicate the magnitudes of nonturbulent atmospheric perturbations that may be of 
concern to sonic-boom propagation studies. Turbulent eddies at the Earth boundary layer distort the 
sonic-boom overpressure time history signature arriving at the ground. Distortions produced by 
turbulent layers far from the ground will tend to return to an undistorted N-wave shape signature 
before reaching the ground. Larger scale wind shears, temperature gradients, and oscillatory activity, 
such as mountain waves and other atmospheric wave motions, may act like airplane maneuvers to 
modulate sonic-boom signatures over distances of a few miles as a result of focusing boom 
overpressure into strong and weak regions. 

Waves are induced by flow over topographic or cloud barriers and their effects are often visible in 
the resulting cloud patterns. Wave phenomena have received much qualitative attention in both 
observation and analysis (refs. 15-17). However, very little statistical documentation exists to 
quantify their intensity, and computational methods have only recently shown the capability to 
realistically handle wave behavior. A particular atmospheric gravity wave model, consisting of a 
solitary wave mode (ref. 18), is proposed here for initial evaluation of nonturbulent atmospheric 
perturbation effects on sonic-boom signatures at ground level. Such a model incorporates realistic 
combinations of vertical and horizontal wind disturbances with their attending temperature or density 
variations. It also approximates the real atmosphere for cases in which a succession of wave crest 
cloud lines are observed downwind of mountains or islands. 


Influence of Atmospheric Variations on Sonic Booms 


• Mean wind and temperature profiles determine 

— Threshold Mach number 
— Boom carpet 
— Intensity effects 

• Small scale turbulence distorts signatures 

• Greater "healing" distances required for quasi-stationary gradients 
and wave motions 

• Assessment of effects on conventional and "shaped" signatures can 
be facilitated by use of discrete models and field measurements 
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REAL-DAY ATMOSPHERIC PROFILE EXAMPLE 


Strong weather fronts and jet-stream activity induce large temperature gradients and wind 
shears in the atmosphere. Horizontal temperature gradients tend to be in balance with vertical wind 
shears, and stable vertical temperature stratification tends to preserve wind-shear layers. Stable 
stratifications occur when the decrease in temperature with altitude (lapse) is less than that 
associated with an adiabatic process following the decrease in pressure with altitude. Thus, 
isothermal layers and inversion layers (where temperature actually increases with altitude) provide 
static stability that tends to maintain wind-shear layers. As the sonic boom propagates through the 
atmosphere, nonuniform changes in the speed of sound along the shock front bend the direction of 
propagation to produce focusing effects. Propagation speeds change directly with wind or air motion 
changes and with temperature. A 1 -percent change in the speed of sound will be experienced for a 
temperature change of approximately 5 °C. 

The figure below illustrates a case (April 10, 1994, Edwards CA) with both strong wind shear in 
a mid-tropospheric stable layer and a strong temperature inversion above the tropopause. 
Windspeeds increase from 23 knots at 13,000-ft altitude to more than 78 knots at 16,000 ft and reach 
peak speeds of 98, 110, and 114 knots at altitudes of 20,000, 25,000, and 34,000 ft, respectively. The 
vector wind shear, given by the wind vector change divided by the altitude interval over which the 
change takes place, reaches a maximum of 0.040/sec at 15,000 ft. Temperature is nearly constant or 
increasing with altitude between 13,300 and 16,500 ft. The vertical temperature gradient or inversion 
rate reaches a maximum of +8 °C/1000 ft between 15,400 and 15,800 ft. 

Above the tropopause, which is marked by a minimum temperature of -66.3 °C at 42,350 ft, the 
temperature increases to -52.3 °C at 45,760 ft. The overall inversion rate for the 12.5 °C temperature 
change in this layer is 3.6 °C/1000 ft. A maximum rate of about 5 °C/1000 ft is observed between 
43,000- and 44,000-ft altitude. The windspeeds begin a marked decrease with altitude above the 
tropopause. 


Edwards Wind and Temperature Profile 

April 10, 1994 




Temperature, deg C 
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WIND SHEAR PROBABILITIES 


Measured wind-shear magnitudes depend on the depth of the layer over which the shear is 
measured. Detailed wind-shear observations and their statistical analyses for altitudes below 
60,000 ft have been abundant since the initiation of the manned Space Program. These observations 
have used special balloons (Jimspheres) tracked by FPS-16 radar to altitudes near 60,000 ft and 
conventional weather balloons to altitudes of about 100,000 ft. Jimsphere wind data studies 
accomplished by the NASA Marshall Space Flight Center (refs. 10-13) have provided a wealth of 
information on wind changes and wind shears and on the relationship between wind shear and the 
thickness of the altitude layer over which the shear is measured. Wind shear is measured to a 
reasonable accuracy over layers thicker than about 300 ft by the Jimsphere and over layers greater 
than about 2000-ft thick by rawinsonde balloons used for routine upper-air weather observations (ref. 
12). Selected data shown below illustrate the variation of the measured wind shear as a function of 
the altitude thickness over which the shear is measured. To bracket a range of stronger shear 
values, the figure illustrates a 90-percentile curve for July, when wind shears are weaker, and a 99- 
percentile curve for February, when they are stronger. Wind shears exceed the values shown by 
these curves for 10 and 1 percent, respectively, of the observations. For convenience of application, 
wind shear in units of speed change per unit altitude distance can simply be expressed in units of 
inverse seconds (or per second). 


Jimsphere Vertical Wand Shear vs Measurement 
Thickness for Two Probabilities 



Wind shear measurement thickness, feet 
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WIND SHEAR PROBABILITIES (cont'd.) 


Wind shears measured over large altitude distances (or thickness scales) have a positive 
correlation with windspeed and therefore tend to decrease from a maximum near jet stream altitudes 
to minimum shear magnitudes in the middle stratosphere, just above supersonic cruise altitudes. 
This decrease in wind-shear magnitudes in the lower stratosphere is depicted below by curves 
extracted from rawinsonde data for altitudes of 45,000 and 80,000 ft (ref. 14). Based on these data 
the maximum shear magnitude measured in the real-day example profile (shown previously) is 
indicated to be beyond the 99.9-percentile value. 


Example Maximum Vertical Wind Shear Probability 
from Rawinsonde Measurements 



Vertical wind shear, per second 
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STRATOSPHERIC TEMPERATURE GRADIENT STATISTICS 


Because the speed of sound in air is a direct function of ambient temperature, the sonic -boom 
propagation direction will bend when strong temperature gradients are traversed at an oblique angle. 
The strongest gradients are associated with inversion layers in the lower stratosphere. Statistical 
analyses of temperature gradient features have not been nearly as adequate as for windspeed and 
wind shear. Reference 14 reports on an examination of a limited amount of rawinsonde data for (a) 
strong lapse rates (temperature decrease with altitude) or (b) strong inversion rates (temperature 
increase with altitude) within layers surrounding the mandatory meteorological reporting levels. This 
study was motivated to improve the knowledge of probabilities for features associated with high- 
altitude turbulence experienced by supersonic cruise aircraft. Values depicting the 1.0-, 0.1-, and 
0.01 -percentile extremes for these data are shown below. 

Note that in the lower stratosphere the magnitudes for the decreasing gradients are not as large 
as those for the increasing temperature gradients. At both the 0.1- and 0.01 -percentile frequency of 
occurrence the decreasing temperature with altitude exceeds the adiabatic lapse rate (3 °C/1000 ft), 
presumably as a result of measurement error in some cases and strong atmospheric dynamics in 
others. The increasing temperature gradient at 0.01 percentile rate, 24 °C/1000 ft, is comparable to a 
density departure rate of approximately 10 percent/1000 ft from the standard day atmosphere, as 
well as to the speed of sound changing by 5 percent/1000 ft from the standard atmosphere values. 


Maximum Temperature Gradients Within Laver 
for Selected Probabilities 


Frequency exceeded 

1.0 % ile 

0.1% ile 

0.01 % ile 

Decreasing temperature gradient, 
°C/1000ft 

-2.9 

-3.6 

-4.2 

Increasing temperature gradient, 
°C/1000 ft 

10 

16 

24 

Temperature increase, 
°C 

8 

10 

12 

Nominal thickness, ft 

for gradients above 5 °C/1000 ft 

1000 

2000 

3000 

Nominal thickness, ft 
for AT above 5.0 °C 

1600 

1000 

500 
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STRATOSPHERIC TEMPERATURE GRADIENT STATISTICS (cont'd.) 


As indicated by the data in the previous figure, temperature increases of more than 10 °C occur at a 
0.1 -percentile rate. Inversion rates of 5 °C/1000 ft extend over layer depths of 3000 ft at a 0.01-percentile 
rate of occurrence. A reciprocal situation, temperature changes of 5 °C (or greater) occurs over layer 
thicknesses of less than 500 ft, also at the 0.01 -percentile rate of occurrence. Probability curves for the 
lapse (decreasing with altitude) and inversion (increasing with altitude) temperature gradients within layers 
of the lower stratosphere are shown in the figure below. 

At times, strong inversion layers are located both above and below strong decreasing temperature 
(lapse) layers in the lower stratosphere. The occurrence of these strong vertical temperature gradients 
exhibits seasonal dependence and indicates some conditional association with high jet stream windspeeds 
and strong lower altitude wind-shear layers. These data indicate that lapse rates in the lower stratosphere 
equal or exceed the dry adiabatic and that inversion rates may exceed 10 °C/1000 ft with overall changes of 
more than 10 °C. In summary, strong inversions can increase the speed of sound by more than 2 percent in 
relatively shallow inversion layers. Depending on wind and flight headings, temperature and wind effects 
may combine to change the speed of sound by nearly 5 percent across layers of 1000 ft and thicker. 


Probability of Maximum Temperature Gradient 

within Laver 
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SINGLE-WAVELENGTH WAVE MODEL 


Atmospheric vertical motion perturbations are subject to wave behavior nearly all the time 
except when the vertical temperature profile does not provide static stability. The resulting wave 
behavior can exhibit multiple modes or wavelengths with both modes and amplitudes changing 
drastically from one altitude layer to another (see outline below). Strong gravity wave motions are 
experienced in both the troposphere and at proposed high-speed civil transport (HSCT) cruise 
altitudes in the stratosphere resulting from dynamic excitation of interactions between the force of 
gravity and atmospheric buoyancy. Resonant interactions are generally possible at all times when 
the temperature decrease with altitude is not adiabatic or superadiabatic. Under favored structures of 
the wind and temperature profiles, barriers to the flow, such as mountains or vigorous cumulus cloud 
lines, and jet stream oscillations can trigger wave motions in which the resulting vertical motions can 
be intense and wave amplitudes can be quite large. Such behavior, with complex upstream 
atmospheric structure, is often beyond solution by classical mathematical analysis. 

Initial analyses of the effects of gravity wave motions on sonic-boom propagation for various 
flight altitudes and trajectories could be overly complicated by the presence of multiple wave modes 
interacting in an atmospheric structure, which changes markedly from one altitude layer to another. 
These complications for initial studies of wave effects on boom propagation and focusing can be 
alleviated by either of two restrictions. 


Attributes of Atmospheric Gravity Wave Behavior 


• Motion perturbations result from Interaction between 

— Wind profile characteristics 
— Presence of flow barriers 

— Temperature profile stability or buoyance effects 

• Wavelength distances between 5 and 50 miles are often observed 

• Large wave amplitudes occur over the full HSCT altitude range 

• Dominant wavelengths generally change with altitude 

• Analytical and numerical treatment of 2-D wave behavior is 

— Fairly extensive and maturing 
— Lacking in field observation validation 

• Realistic, simplified representations of wave behavior in the 
troposphere and lower stratosphere can be generated 
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SINGLE-WAVELENGTH WAVE MODEL (cont'd.) 


One method is to restrict attention to two-dimensional wave cases that produce high-amplitude 
perturbations for only one or one-and-a-half wavelengths. Another means of simplification is to 
restrict attention to two-dimensional wave cases that produce a single-mode (i.e., single- 
wavelength) trapped wave. This second simplification is believed of more general interest and is 
described below. 

The schematic below shows a two-dimensional trapped, single-mode wave to illustrate the 
variation of wave induced motion perturbations in the various wave sectors. The wave is produced 
by a numerical simulation (ref. 18 ) using an analytically selected upstream profile with constant 
shear (depicted on the left) and stability (i.e., constant lapse rate). Flow following the streamlines 
experiences the greatest vertical motion perturbations at the altitude of maximum wave amplitude 
and at inflection points between troughs and crests. Vectors (not to scale) show downward 
perturbation velocities from the crest to the trough and upward velocities from the trough to the 
crest. Temperature perturbations also reach their greatest amplitudes at the altitude of maximum 
wave amplitude. Temperature increases (warming) occur in the troughs where the air has descended 
to higher pressures with adiabatic heating; temperature decreases (cooling) occur in the crests 
where the air has ascended from its upstream altitude level. Horizontal wind perturbation vectors 
(not to scale) show that the wave dynamics increase horizontal windspeeds in the lower section of 
the troughs and the upper sections of the crests. Horizontal speeds are reduced in the upper section 
of the troughs, where the wave amplitude decreases with altitude, and in the lower section of the 
crests, where rotors with reverse flow near the ground are sometimes created. 


Gravity Wave Schematic with Wave-Induced 
Perturbation Velocity Components 


Wind profile 


Altitude 



Streamline field 
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MODEL WAVE INTENSITIES 


Selected barrier conditions and atmospheric structure were used to produce a single-mode wave, as 
described above, by both analytical solution and numerical simulation (ref. 18, Case II). This case 
generated a nominal intensity wave using a barrier height of approximately 1600 ft in a wind of 20 knots at 
die base altitude and increasing a modest 2 knots/1000-ft altitude. The temperature lapse rate was 2.25 
°C/1000 ft, fairly typical of the lower troposphere and near the standard atmosphere value of about 2 
°C/1000 ft. Maximum wave amplitude is attained at an altitude near 10,000 ft, where the stream function's 
trough-to-crest height change is approximately 4600 ft and the temperature change from trough to crest is 
about 3 °C. At 10,000-ft altitude the upstream windspeed is about 40 knots, and a wavelength of about 11 
nmi prevails throughout the wave pattern. Horizontal windspeed varies by 16 knots from maximum to 
minimum at the surface and 3 to 4 knots at 15,000-ft altitude. The maximum up- and downdraft velocities 
are approximately 8 ft/sec. Strength of this wave would be subjectively rated as light or weak, on the basis 
of the induced vertical velocity and its limited extent above the barrier height Its amplitude and wavelength 
are fairly representative of extended (or trapped) wave conditions in which several lines of crest clouds are 
observed downstream of flow barriers. A stronger wave case, in which the maximum vertical velocity is 
greater by a factor of at least 3, should also be used in boom propagation studies to encompass a large 
portion of likely wave amplitudes. 

Strong temperature gradients and wind shears are generated locally within the wave structure. The 
boom propagation effects of similar strength shears and gradients in quasi-continuous layers should be 
compared with the wave cases to delineate the added focusing effects of vertical motion and wave shape. 
Layer slopes with respect to flightpath could be varied to cover the range experienced in the wave 
structures. Based on climatology and clear air turbulence experience, intensities of wind shear used for 
boom propagation studies should range from less than 0.010/sec (light) to more than 0.020/sec (strong). A 
realistic range for temperature gradients, which could combine with the wind shear, is from -3 °C/10Q0 ft 
to more than +10 °C/1000 ft. 


An Approach for the Study of Boom Propagation 
through Atmospheric Waves 

• Lower- and mid-tropospheric wave cases 

— Two representative wave intensities or amplitudes 
— Flight altitudes from level of max amplitude to cruise 
— Mach range at each altitude to cover most aircraft types 

• Upper-troposphere wave cases 

— Two wave intensities and flight conditions varied as above 
— Limit altitude band of wave excitation 
— Use both up- and down-wind headings 

• Lower-stratosphere wave cases 

— More cases to cover greater range of wind and temperature profiles 
— Two wave intensities and flight conditions varied as above 
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MODEL WAVE INTENSITIES (cont'd.) 


Variables in the boom-propagation study will include flight altitude and Mach number as well 
as the wave models. Suggestions for the scope of conditions to be studied are listed below. The 
generation of simple, yet realistic, wave models in the upper troposphere and lower stratosphere will 
require added consideration of the input barrier and atmospheric structure. 


CLOSING REMARKS 


Atmospheric gradients and wave motions can affect sonic-boom propagation. Simple, discrete 
examples that can guide studies of boom propagation have been described in this paper. The 
specification of wave cases for boom work should initially avoid the potential complexities that could 
result from variation of wave amplitude and wavelength from one altitude layer to another. It is 
recommended that initial studies of boom propagation through wave perturbations be attended by 
propagation cases through similar-strength stratified layers that slope in relation to the flightpath. 


• Sonic boom propagation is subject to atmospheric turbulence, 
gradients and waves 

• Complex atmospheric wave behavior can easily complicate 
boom propagation studies 

• Realistic simplifications of wave behavior will facilitate the 
interpretation of boom propagation study resutis 

• Comparison of wave induced boom focusing effects with those 
from stratified gradients of similar intensity is recommended 
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USAF SINGLE-EVENT SONIC BOOM 
PREDICTION MODEL: PCBoomS 


ABSTRACT 

The Air Force has developed PCBoom3, a general-purpose, single-event sonic boom prediction 
model. The model operates on an IBM PC or compatible, under DOS or Windows. It is accessed via an 
integrated environment which controls building of input cases, running boom calculations, displaying 
contours and signatures, and managing all associated data. The primary boom calculation is via a variation 
of FOBOOM, the focus-boom extension of Thomas's program. Arcraft input is either via a user-provided 
F-fiinction, or simple N-wave F-functions tabulated for about 20 current aircraft. A fast boom calculation, 
based on Plotkin's SBORT algorithms, is included for simple N-wave F-functions in a windless atmosphere 
and flight altitudes up to 60,000 feet. After a run is complete, the user can access an index identifying 
significant events (focal zones, beginning of footprint, etc.), then plot boom amplitude contours and 
signatures or spectra at any point in the footprint. The primary uses of this program are expected to be 
operational planning and boom incident investigation. However, because of the commonality between 
FOBOOM and the MDBOOM program currently being used for low boom configuration design, this 
program is of interest to the HSCT community, especially as supersonic route planning activity increases. 

The Ar Force recently conducted a flight test program to evaluate the focal zone capabilities of 
PCBoom3 . Initial results of that program validate the prediction of focal zone geometry, amplitudes, and 
waveforms. 

[This work was sponsored by USAF AL/OEBN.] 
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PCBoom was developed to satisfy three basic Air Force needs. The first is a tool to predict sonic 
booms for use in environmental assessments of proposed actions. The second is a planning tool to 
minimize impact to sensitive areas. An example of this type of planning is if the operators of an offshore 
supersonic range would like to set daily minima for how far out the range users must go. The third 
application is to quantify what happened when an unintended incident does happen. 

These needs lead to the requirement for a program that is relatively easy to use, runs quickly, and 
can be hosted on a PC. It is essential that input is as flexible as possible, and that output be available in 
graphical format. 


AIR FORCE REQUIREMENTS 

• Single-Event Sonic Boom Impact Prediction 

• Planning Tool 

• Sonic Boom Incident Investigation 
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The program has several methods for specifying trajectories. First, there is a file structure similar 
to that of FOBOOM/MDBOOM. Second, there is a general "Maneuver Driver" which allows interactive 
description of a maneuver. This is aimed at the planner or airspace manager who must be able to work 
with virtually any description of a maneuver. The interactive form of "ManD" helps the user reconcile 
information that may be inconsistent. Third, the program can import trajectory data from sources such as 
radar tracking, ACMI, or pre-computed trajectories. The program is designed to import these from a 
relatively benign ASCII file format. A user with unique data forms can prepare a conversion routine to this 
specification. 

The program uses standard sonic boom ray tracing theory, and includes focus analysis via Guiraud's 
similitude and the Gill-Seebass numeric focus signature. Use of this program for environmental planning 
and claim investigation requires that it employ state-of-the-art methods. 

Outputs are contours of equal overpressure (or other metrics), isopemp charts which provide a 
good visual interpretation of footprints and focal zones, and signatures and spectra. 

GENERAL FEATURES 

• General Flight Profile Input 

• Full Ray Tracing Scheme 

• Non-Standard Atmosphere, With Winds 

• Focus and Post-Focus Boom Impact Regions 

• Outputs: Contours, Footprints (Isopemps), Signatures, 

Spectra 
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The user works primarily with a simple pull-down menu system. The goal of this module is to set 
up the data stream needed by the computational modules that actually do the work. When a computational 
module is ready, MIlVl writes a batch file to execute it, then exits. The batch file ends by re-running MIM, 
putting the user back in the menus. This program is PC specific, but would be fairly simple to replicate on 
other interactive systems. 


PROGRAM STRUCTURE: 

MASTER INPUT MODULE (MIM) 

• Manages All Data, User Inputs, Via Simple Menu System 

• Sets Up Runs Stream for Computational Modules 

• Initiates Computational Modules as Separate Processes 

• PC-Specific Interface 



The real work is done in the computational modules. These are written in Fortran 77, and for the 
most part are portable. Graphics are done with a PC-specific Fortran-callable plotting library with 
"industry standard" call styles. ManD, PCBPlot, and SIGOUT have PC-specific interactive interfaces, but 
these are kept isolated from the computational parts. 

Boom calculations are done by either the full-theory FOBoom3 module or, for applications 
requiring very fast analysis, a scheme based on the rapid SBORT algorithms. Both schemes yield a file of 
boom signatures on the ground. This file is processed by a "footprint processor" which organizes it and 
indexes it for use by the output modules PCBPlot and SIGOUT. 

It is worth noting that the FOBOOM3 module is derived from the same code from which 
MDBOOM was developed. It does not contain the proprietary configuration analysis tools, but those have 
been removed by blocking out subroutine calls and removing the corresponding modules. MDBOOM 
users may find it useful to adapt some of the new features to their work. As many routines as possible 
have been kept identical between the two programs, which will be of long-term benefit to maintenance 
of both. 


COMPUTATIONAL MODULES 

• Maneuver Driver (ManD) 

• FOBoom3 Boom Calculator 

• SBORT Simple Boom Calculator 

• Footprint Processor 

• PCBPlot Footprint/Contour Output 

• SIGOUT Signature and Spectrum Output 

• All Written in Fortran 77 

• "Industry Standard" Plotting Calls 
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The contour plotting module will draw isopemps and contours. The isopemps - the ground 
intercepts of boom generated at specific times - are very useful for seeing where the boom footprint lies, 
and in interpreting focal zone geometry. The contours quantify the boom amplitude. Loudness is 
computed by the methods used by NASA-Langley, with a 125 msec auditory time constant. 


PCBCont CONTOUR/FOOTPRINT MODULE 

• Isopemps (Ray/Ground Intersections) 

• Contours of P max , CSEL, ASEL, PL 

• Interactive User Control of Plot Format 
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This is an isopemp plot for a supersonic turn at constant Mach number. There is a "+" marking the 
beginning of the calculated maneuver, and "o"s marking the entry point. The caustic location during the 
turn is apparent, and the multi-valued region around the turn entry cusp can be seen. 


TYPICAL FOOTPRINT ("ISOPEMP") OUTPUT 
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This is the corresponding overpressure contour. 


TYPICAL CONTOUR OUTPUT 
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The signature output module allows plotting signatures and spectra at an arbitrary location in the 
footprint. It interpolates signatures between the actual computed grid points. It is aware of multiple sheets 
in focal regions, and the absolute phase, so it can assemble complex N-U signatures seen near focal zones. 
The module plots energy density spectra and residual shock spectra. Spectra can be either narrow or one- 
third octave band. The one-third octave band spectra are normalized by the 125 msec auditory 
time constant. 


SIGOUT SIGNATURE/SPECTRUM MODULE 

• Plots Time History Signatures 

• Complex M N—U" Signatures in Focal Zones 

• Spectra and Residual Shock Spectra 


* 
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This is an example signature output. This particular one is the maximum focus condition under the 
flight track for an F-16 performing maximum-power level acceleration at 10,000 feet. The peak pressure is 
17 psf, and the CSEL is 121.5 dB. 


Calculated Focus Signature 
F-16 Level Acceleration at 10,000 -feet MSL 



Time, milliseconds 




Focal zone flight tests were recently conducted. An array of BEARS with 500- to 1,000-foot 
spacing was placed across calculated focal zones. Observers were also located along the arrays, and the 
audible sensations of "distant thunder" on the shadow side, "BANG-BANG" at the focus, and "carpet 
boom with afterpops" filled expectations. There were 49 passes over the arrays, with a good mix of 
maneuver types. Except for a very few anomalous runs that missed altogether, all of the foci intercepted 
the expected locations. The behavior away from the focal zone itself- with N and U waves separating - 
agreed well with predictions. Boom amplitudes also agreed with predictions. 


FLIGHT TEST VALIDATION: 

Project Have BEARS 

. USAF AL/OEBN and TPS effort, April 1994, Edwards AFB 

• F-16 Performing Level Acceleration, Diving Acceleration, Steady Turn, 
Pullup/Pushover Focusing Maneuvers 

• Signatures Recorded With Arrays of USAF BEARs 

• Excellent Results for Focal Zone Location, Geometry, Boom Amplitude 
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This is a typical measured boom at the focus condition. This is for a full-throttle acceleration, 
corresponding to the calculated signature shown earlier. The peak pressure is 17 psf, and the CSEL is 
123.8 dB. Both are in very good agreement with prediction. The shape differs somewhat, and there are 
spectral differences, so the scaling of the Gill-Seebass solution needs to be examined, but the result is 
generally quite good. Having the recorded data in digital form from the BEARs provides the first real 
opportunity for a detailed comparison of theory with flight test. 


File B19124QR.413 12:40:32.61 April 13 1994 

Pm3x = 17.08 Pmin = -4.54 16001 points Site 19 S/N 4019 



0 . 50 . 100 . 150 . 200 . 250 . 300 . 350 . 

Time, milliseconds 
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CONCLUSIONS 


PCBoomS Is Nearing Completion 

Provides Simple Access to State-of-the-Art Sonic Boom 
Prediction 

Validation by Recent Flight Tests 

Commonality With MDBOOM Will Benefit the HSCT 
Community 

Expected Release in Fall 1994. Will be Available From 
AL/OEBN 
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INTRODUCTION 


The inability to systematically predict community response to exposure to sonic booms (and other 
high energy impulsive sounds) is a major impediment to credible analyses of the environmental effects 
of supersonic flight operations. Efforts to assess community response to high energy impulsive sounds 
are limited in at least two important ways. First, a paucity of appropriate empirical data makes it 
difficult to infer a dosage-response relationship by means similar to those used in the case of general 
transportation noise. Second, it is unclear how well the "equal energy hypothesis" (the notion that 
duration, number, and level of individual events are directly interchangeable determinants of 
annoyance) applies to some forms of impulsive noise exposure. 

Some of the issues currently under consideration by a CHABA working group addressing these 
problems are discussed. These include means for applying information gained in controlled exposure 
studies about different rates of growth of annoyance with impulsive and non-impulsive sound exposure 
levels, and strategies for developing a dosage-response relationship in a data-poor area. 
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The state of the art in prediction of community response to high energy impulsive noise exposure - 
essentially, sonic booms, artillery fire, and blasting - is much less advanced than prediction of 
community response to general transportation noise. The most obvious difference between the two 
cases is between an embarrassment of social survey information about the annoyance of general 
transportation noise and a scarcity of similar findings about the annoyance of impulsive noise. 

Fields (1991) notes that more than 300 social surveys of community response to the general din of 
the urban environment have been conducted in the last few decades. Schultz (1978) found enough 
quantitative information about 161 of these data points to permit his well-known synthesis of a 
descriptive dosage-response relationship between Day-Night Average Sound Level and the prevalence 
of a consequential degree of self-reported annoyance in urban settings. Limiting themselves to the 
same selection criteria as Schultz (1978), Fidell, Barber and Schultz analyzed another 232 data points 
in 1991. Others (notably FICON, 1992) have suggested alternate fitting functions for various subsets 
of the same data, as summarized in Table 1 of Fidell and Pearsons (1993). 

As seen in Figure 1, information about community response to general transportation noise is now 
so abundant - and so variable - that it is doubtful that information collected in any further social 
surveys can substantially affect commonly accepted fitting functions and interpretations. As a 
complement to purely descriptive accounts of this mass of data, efforts such as that of Green and 
Fidell (1991) have recently begun to develop theory-based explanations for the variability in annoyance 
rates in different communities with similar noise exposure. 



Figure 1 Social survey findings on the prevalence of 
annoyance associated with residential exposure to general 
transportation noise. (Fitting function from one parameter 
model of Green and Fidell, 1991.) 


186 



In contrast, only a few social surveys have been conducted on the effects of high energy impulsive 
noise exposure on communities, and only one well-accepted synthesis of a descriptive dosage-response 
relationship for these data has been completed. Figure 2 summarizes the information presently 
available from which a dosage-response relationship between impulsive noise exposure and the 
prevalence of annoyance can be directly synthesized. The data displayed in Figure 2 are derived from 
the well-known Oklahoma City study (Borsky, 1966); from two surveys of reactions to artillery noise 
(Schomer, 1985); and from a recent NASA study of reactions to sonic booms reported elsewhere in 
these proceedings by Fields, Moulton, Baumgartner and Imm-Thomas. All of the original noise 
measurements at Ft. Lewis and Ft. Bragg were modified to account for pressure doubling. The plotted 
values also reflect new information about the Oklahoma City noise measurements (Schomer, 1994). 

The best known dosage-response relationship (that of CHABA Working Group 84) was not based 
on even this much information, however. Galloway (1981) was forced to develop a dosage-response 
relationship from only 14 data points. Another CHABA Working Group on Assessment of 
Community Response to High Energy Impulsive Sounds now revisiting the problem 13 years later is 
attempting to interpret what has been learned empirically and theoretically about community response 
to impulsive noise exposure in the interim. 

The most obvious remedy for the paucity of information about community response to impulsive 
noise exposure is the conduct of several new large scale social surveys. This is not a realistic 
possibility for several reasons. First, opportunities to conduct surveys of adventitious exposure to high 
energy impulsive noise are limited. Because no civil aircraft fly supersonically over land, and because 
the military confines its overland supersonic operations to areas of relatively low population density, 
only small populations have yet experienced high levels of exposure to sonic booms on a regular, long 



Figure 2 Social survey findings on prevalence of 
annoyance associated with exposure to high energy 
impulsive noise. 
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term basis. Blasting, artillery training and other sources of non-aircraft, high energy impulsive noise 
are also generally experienced only by localized populations. 

Second, opportunities for conducting field studies involving intentional exposure of entire 
communities to sonic booms have been all but foreclosed in the United States since passage of the 
National Environmental Policy Act of 1969. 

Third, it is doubtful that information about effects of sonic booms in the exposure ranges of 
greatest interest for purposes of developing a dosage-response relationship can ever be collected. C- 
weighted DNL values in systematic studies of community response to high energy impulsive noise 
exposure that have been conducted to date are all in the region below 70 dB. No community has yet 
experienced the numbers of daily supersonic flights necessary to produce greater long term impulsive 
noise exposure, nor is it possible for a variety of reasons to create or credibly simulate such exposure 
on a large scale. 

Controlled studies of subjective judgments of the annoyance of individual impulsive sounds can 
contribute information that might be used to complement social survey findings. Kryter, Johnson and 
Young's (1968) field study comparing the annoyance of sonic booms to that of subsonic aircraft 
overflights is one source of such information that was considered by CHABA Working Group 1984. 
Schomer (1994) has since collected newer information about the relative annoyance of other impulsive 
sounds. This information was derived from direct paired comparison judgments of the annoyance of 
impulsive and non-impulsive noises heard in laboratory and field settings. Schomer’s outdoor noise 
measurements of sounds presented for judgment in field testing were made at the time of signal 
presentations, and pressure doubled for the sake of consistency with those typically made of sonic 
booms. 

In Schomer's data set, the annoyance created by an impulsive noise event of a given CSEL grows 
at twice the rate as the annoyance of a non-impulsive noise event of a numerically equivalent ASEL. 
This means that CDNL values cannot be directly calculated for purposes of predicting annoyance as a 
logarithmic sum of CSEL values of constituent noise events. Just as OSHA employs a 5 dB/doubling 
rale rather than a 3 dB/doubling rale for calculating noise exposure for purposes of predicting hearing 
damage risk, Schomer' data suggest that a 6 dB/doubling rale is required when CDNL values are used 
for purposes of predicting the prevalence of annoyance due to a summation of high energy impulses. 

Alternatively, a greater rate of growth of annoyance for impulsive than non-impulsive noise events 
could be interpreted as requiring that the slope of dosage-response relationship for community response 
to the impulse noises be considerably steeper than the slope of the dosage-response relationship for 
non-impulsive noise. 

One might think that a difference in the fit to the data of two dosage-response relationships with 
slopes differing by a factor of two would be immediately apparent from simple visual inspection. This 
is not the case, however, because the small number of data points and their considerable variability do 
not greatly constrain the shape of a fitting function. It is therefore important to explore other means 
for checking the reasonableness of drawing inferences for the shape of a dosage-response relationship 
for community response to impulsive noise exposure from the findings of controlled exposure studies. 
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One reasonableness check can be made by analogy with the case of predicting community response 
to non-impulsive noise. Green and Fidell (1991) attempted to infer the shape of a dosage-response 
relationship from first principles, rather than through statistical curve fitting exercises. They likened 
the prevalence of annoyance in a community to a response to a dose of noise exposure: like any other 
treatment administered to a human population, the response to the same dose can be expected to vary 
from individual to individual, and hence from community to community. 

Green and Fidell suggested that the relationship between dose and response has two components: 
one associated with the effective loudness of the noise exposure, and one associated with the sum of 
all nonacoustic influences on self-reports of annoyance. The nonacoustic influences may be considered 
in the aggregate as a form of response bias. The former term establishes the slope of the dosage- 
response relationship, whereas the latter establishes its position along the abscissa. 

Green and Fidell showed that a one-parameter model provides a good account of the relationship 
between exposure to general transportation noise and the prevalence of annoyance in communities. It 
is thus reasonable to ask whether the same model can also be applied to the case of impulsive noise 
exposure. The one parameter model assumes that reactions of community members to noise exposure 
are exponentially distributed with a mean population value, m. The value of m is assumed to be 
related to the Day-Night Average Sound Level by: 

10 log m = 0.3 L dn 

Thus, noise exposure creates a distribution of reactions within a community with a mean value that 
increases with the level of noise exposure. Individuals describe themselves as highly annoyed when 
their reactions to noise exposure exceed a criterion value for reporting annoyance. The proportion of 
the population describing itself as highly annoyed is predicted as 

p __ ^-(A/m) 

where P is the probability of reporting high annoyance, m is defined as above, and A is the criterion 
value for reporting annoyance. 

Figure 3 shows the fit of the Green and Fidell model to the data displayed in Figure 2. The curve 
is generated by the relationship shown in Equation 1. The value of the parameter that controls the 
slope of the curve - that is, the power to which DNL is raised to calculate the effective loudness of 
noise exposure - is 0.3. The horizontal position of the curve on the abscissa is determined by the 
average value of CDNL at which respondents describe themselves as highly annoyed in this data set. 
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Figure 3 Fit of one parameter model (Green and Fidell, 
1991) to data displayed in Figure 2, exponent = 0.3. 


A predictive relationship that provides a useful account for a data set should reduce the amount of 
variance unaccounted for in the data set. Figure 4 shows that this is in fact the case by replotting the 
data of Figure 3 in a manner that removes the effects of response bias. The abscissa of Figure 4 
subtracts the average response bias observed in each study (expressed in units of CDNL) from the 
measured CDNL values for each data point. This normalization of each data set to the average value 
of CDNL at which survey respondents in a study described themselves as highly annoyed effectively 
removes response bias, and leaves apparent the correspondence between the observed and predicted 
rates of growth of annoyance with exposure level. The reduction in variability in Figure 4 with respect 
to that of Figure 3 is readily apparent. 

If annoyance with impulsive noise exposure grows with level at a rate faster than it grows for non- 
impulsive noise, one might expect a larger exponent for CDNL than for ADNL. Schomer (1989), for 
example, has suggested that a value of 0.4 might be more appropriate than 0.3. Figure 5 shows the fit 
of the Green and Fidell model with a value of 0.4 for its one free parameter to the same data set. 

The standard deviation of the differences between the observed and predicted proportions highly 
annoyed for the fit shown in Figure 3 is 3.2%, while the corresponding standard deviation for the fit 
shown in Figure 5 is 3.8%. The differences in standard deviations are so slight that there is no 
compelling argument for adopting one that suggests a higher rate of growth of annoyance with 
exposure level in the impulsive case. 

Thus, it is not yet apparent how the observations of controlled exposure studies in which the 
annoyance of impulsive and non-impulsive sounds are directly compared can be applied to derivation 
of a dosage-response relationship for the prevalence of annoyance due to high energy impulsive 
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sounds in communities. Further analyses of these data are expected before the current CHABA 
working group makes any new recommendations. 



| A OK CITY X NASA □ FT. LEWIS m FT BRAGG | 

Figure 4 Reduction in variability in social survey findings 
displayed in Figure 3 attainable by normalization to remove 
effects of response bias. 



Figure 5 Fit of one parameter model (Green and Fidell, 
1991) to data displayed in Figure 2, exponent = 0.4. 
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INTRODUCTION 


The proposed development of a second-generation supersonic commercial transport has resulted in 
increased research efforts to provide an environmentally acceptable aircraft. One of the enviromental 
issues is the impact of sonic booms on people. Aircraft designers are attempting to design the transport to 
produce sonic boom signatures that will have minimum impact on the public. Current supersonic 
commercial aircraft produce an "N-wave" sonic boom pressure signature that is considered unacceptable 
by the public. This has resulted in first-generation supersonic transports being banned from flying 
supersonically over land in the United States, a severe economic constraint. By tailoring aircraft volume 
and lift distributions, designers hope to produce sonic boom signatures having specific shapes other than 
"N-wave" that may be more acceptable to the public and could possibly permit overland supersonic flight. 
As part of the effort to develop a second-generation supersonic commercial transport, Langley Research 
Center is conducting research to study people's subjective response to sonic booms. As part of that 
research, a system was developed for performing studies of the subjective response of people to the 
occurrence of simulated sonic booms in their homes. 

The In-Home Noise Generation/Response System (IHONORS) provides a degree of situational 
realism not available in the laboratory and a degree of control over the noise exposure not found in 
community surveys. The computer-controlled audio system generates the simulated sonic booms, 
measures the noise levels, and records the subjects' ratings and can be placed and operated in individuals' 
homes for extended periods of time. The system was used to conduct an in-home study of subjective 
response to simulated sonic booms. The primary objective of the study was to determine the effect on 
annoyance of the number of sonic boom occurrences in a realistic environment. 
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IN-HOME NOISE GENERATION/RESPONSE SYSTEM DIAGRAM 


A diagram of the IHONORS system is shown in figure 1. The system consisted of a computer and 
compact disc player that played the simulated sonic booms at randomly-selected, pre-programmed times 
through a preamplifier and amplifier into three or four loudspeakers located in different rooms of the 
house. The two indoor microphones and sound level meters measured the levels of the booms as they 
occurred and also continuously measured the ambient noise levels in the home. The measurements were 
then transferred to the computer and stored on its hard disk. At the end of the day, the test subject used 
the trackball to answer a series of questions about his or her activities during the day and his or her 
subjective response to the sonic booms heard. Once a week the noise measurements and test subject 
responses were downloaded via a modem to a central computer. The data were then checked to ensure 
that the system was operating correctly. 
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Figure 1. 





DEPLOYED IN-HOME NOISE GENERATION/RESPONSE SYSTEM 


Figure 2 shows the actual components as deployed in a test subject's home. Selected homes were 
limited to single family detached dwellings so that the simulated sonic booms would not penetrate the walls 
into neighbors homes as might be the case in apartments or townhouses. The range of homes selected 

covered the economic range from lower middle class to upper middle class. The computer, compact disc 
player, trackball, and monitor were placed on the shelves of a microwave oven cart. The preamplifier, 
amplifiers, and sound level meters were placed inside the lower cabinet of the cart. The cart was placed in 
a position that was convenient for the test subject and that simplified the installation of cables. Three or 
four loudspeakers were placed in the rooms that the test subject indicated he or she most used during the 
14 hour boom period each day. The two microphones were placed in two of the rooms with loudspeakers. 
The system components blended in well with the existing decor in most cases and were often decorated 
with bric-a-brac by the test subject. 


2 Indoor Microphones Test Subject's Home 3-4 Loudspeakers 



Computer, Track Ball. Modem, CD Player, Amplifiers, Sound Level Meters 


Figure 2. 
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TEST PLAN 


The final test plan is outlined in figure 3. Eight IHONORS systems were used to conduct the in- 
home study of subjective response to simulated sonic booms. A system was deployed for eight weeks in 
each of 33 homes. Each day the system played simulated sonic booms during a 14 hour period as the test 
subject went about his or her normal activities. At the end of the 14 hours the test subject rated his or her 
annoyance to the sonic booms heard during the day. A total of 264 weeks of data including over 1800 
subjective annoyance responses to daily sonic boom exposures were collected, the equivalent of five years 
of sonic boom exposures in realistic environments. 


. February 1 993 to December 1 993 
* 33 homes 


• 14 hour test day (no booms during normal sleep period) 


1848 total exposure days 


58,443 total sonic booms 


Figure 3. 
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EXPERIMENT DESIGN 


As shown in figure 4, the sonic booms presented each day represented combinations of three sonic 
boom pressure signatures or waveforms, three A-weighted sound exposure levels (SEL(A)), and seven 
sonic boom frequencies. The SEL(A) levels used were nominally 66, 70, and 74 dB. These values 
covered the range of indoor sonic boom levels estimated for a variety of second-generation supersonic 
transport designs. The pressure waveforms represented an outdoor N-wave, an indoor N-wave, and an 
outdoor "shaped" wave. All the sonic booms had a rise time (x) of four msec and a duration of 300 msec. 
The frequencies were 4, 10, 13, 25, 33, 44, and 63 booms per 14 hour period. Only one sonic boom 
waveform was presented each day. On most days the sonic boom was presented at only one SEL(A) level. 
On a few days the sonic boom was presented at two or three of the SEL(A) levels. 


. 3 pressure signatures 

. outdoor N-wave, T = 4 msec 

. indoor N-wave, T = 4 msec 

. outdoor shaped, T = 4 msec 

. 3 levels - 66, 70, 74 dB - SEL(A) 

. 7 boom occurrence rates - 4, 10, 13, 25, 33 ,44, 63 booms 

per day 

Figure 4. 
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TEST SUBJECT RESPONSES 


The information obtained from the test subjects is outlined in figure 5. The computer-generated 
questions answered by the test subject at the end of each day are summarized as follows: (1) when were 
you not inside your house, (2) what activities did you do while inside your house, (3) how annoying were 
the sonic booms you heard today on a 0 to 10 scale, and (4) were you startled by any of the sonic booms 
today? In addition to the daily questions, each test subject answered pre- and post-test questionnaires 
similar to those used in studies surveying people exposed to aircraft-generated sonic booms. 


. Pre- and post-test questionnaires 

. Daily computer-generated questions 

. When were you not inside the house? 

. What activities did you do while in the house? 

. How annoying were the sonic booms you heard 
today? (0 to 1 0) 

. Were you startled by any of the sonic booms today? 
(Yes or no) 


Figure 5. 
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TEST DAY SUMMARY 


The eight IHONORS systems were deployed in 33 test subjects' homes for a total of 1848 exposure 
days. Figure 6 gives a breakdown of how the test subjects responded during those 1848 days. Over 90 
percent of the time, the subjects responded with an annoyance judgment. This was better than expected. 
Subjects were outside of the home for the entire 14 hour period only one percent of the time. Subjects 
were out of town or forgot to respond to the questions less than three percent of the time. Technical 
problems with the system caused the loss of data on less than six percent of the days. Although the 
systems were designed to and did recover by themselves from power outages by rebooting, the most 
common technical problem was a momentary brown-out or flicker in the power supply causing the 
computer to lose its default memory settings but not triggering a reboot of the system. 


Type of day 

Number 

Percentaae of total 

Subject responded with annoyance judgment 

1676 

90.7 % 

Subject was outside all day 

19 

1 .0 % 

Subject did not respond 

47 

2.6 % 

Day voided due to technical problem 

106 

5.7 % 


Total exposure days 


1848 


100 % 


Figure 6. 
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EFFECT OF NUMBER OF SONIC BOOM OCCURRENCES ON ANNOYANCE 


Figures 7 and 8 illustrate the preliminary results concerning the effect of number of booms heard. 
Figure 7 shows subjective annoyance ratings versus the SEL(A) level of the individual sonic boom repeated 
during the day. The number of occurrences (n) of the boom is divided into five intervals, each having 
roughly the same number of data points. The linear regression lines for each interval are plotted in the 
figure. As illustrated in the figure, the subjective annoyance response increases as the number of 
occurrences of a sonic boom increases. 


Subjective 

annoyance 

rating 



Figure 7. 
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MODELING THE EFFECT OF NUMBER OF SONIC BOOM OCCURRENCES ON ANNOYANCE 


This effect of number of occurrences can be modeled by the addition to the measured sonic boom 
level of the term "k * log(number of occurrences)" : Figure 8 shows the values of "k" and the 
corresponding 95% confidence intervals determined from regression analyses for each of several noise 
metrics. The metrics were perceived level (PL), two versions of Zwicker's loudness level (LL^d and 
LLz f)> perceived noise level (PNL), A-weighted sound exposure level (SEL(A)), C-weighted sound 
exposure level (SEL(C)), and unweighted sound exposure level (SEL(U)). The metrics represent different 
ways of modeling the frequency response of the ear (refs. 1,2). As shown in figure 8, analyses of the data 
indicate that the calculated value of "k" ranged from 10 to 15 depending on the noise metric considered. 
However, for almost all the metrics, the 95% confidence interval about the calculated value includes the 
value of 10. Therefore, in those cases, the value of "k" cannot, statistically, be said to be significantly 
different from 10, the predicted value based on energy addition. 

The two metrics whose 95% confidence intervals on "k" did not include the value of 10 were 
examined further. The value of "k" was calculated from the regression coefficients aj and for each 
metric. Comparison of a2, the "log(number of occurrences)" coefficient, across all metrics found that the 
value of a2 was constant across all metrics, indicating that the differences in "k" for SEL(U) and SEL(C) 
were due to the metric and not the effect of number of occurrences. 

Annoyance = ao + a-| * Level + a2 * log(Number of occurrences) 



Noise Metric 

Figure 8. 
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COMPARISON OF NOISE METRICS 


Having confirmed the model for summing the effect of multiple sonic booms, the total daily sonic 
boom exposure was calculated in terms of each of the noise metrics for comparison with the test subjects' 
daily annoyance judgments. Daily exposure is commonly expressed in terms of Day-Night Level (DNL). 
Day-Night Level uses "10*log(number of occurrences)" to sum multiple events and then averages the noise 
energy across 24 hours (ref. 2). (No late night penalties were assessed.) Although DNL is usually 
associated with A-weighted sound pressure level, for the purpose of this study, a DNL was calculated 
using each of the noise metrics. 

The predictive ability of the noise metrics was then compared based on the correlation coefficients 
between the subjective annoyance ratings and the different DNL's. The resulting rank order of metrics is 
shown in figure 9. Perceived level was the best predictor of annoyance to the simulated sonic booms. 
Figure 9 also illustrates the subjective annoyance ratings as a function of DNL(PL) and as a function of 
the more commonly used DNL(A) and DNL(C). A best-fit regression line is drawn through the data in 
each plot. 



Figure 9. 
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EFFECT OF SONIC BOOM WAVEFORM SHAPE ON ANNOYANCE 


Figure 10 shows the effect of sonic boom waveform on annoyance for noise measurements based 
on perceived level. The figure shows the regression line of subjective annoyance rating on day-night level 
based on perceived level for each of the three sonic boom waveforms used in the study. Indicator 
(dummy) variable analyses found no difference in either slope or intercept between the three different sonic 
boom waveforms. This result was the same for all the noise metrics considered. 


Subjective 

annoyance 

rating 



Figure 10. 
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EFFECT OF STARTLE ON ANNOYANCE TO SONIC BOOMS 


Figure 1 1 shows the effect of startle on annoyance for noise measurements based on perceived 
level. The regression lines for both startled and not startled are shown for subjective annoyance rating 
plotted against day-night level based on perceived level. Indicator (dummy) variable analyses indicated a 
significant difference in both slope and intercept between startled and not startled by a sonic boom. 
Annoyance is greater when the test subject is startled and the magnitude of the increase in annoyance 
increases as the sonic boom exposure increases. This effect was the same for all the noise metrics 
considered. 


Subjective 

annoyance 

rating 



DNL (PL), dB 


Figure 11. 
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PERCENTAGE OF TEST SUBJECTS HIGHLY ANNOYED AS A FUNCTION OF A-WEIGHTED 

DAY-NIGHT LEVEL 


The impact of aircraft flyover noise has usually been examined in terms of the percent of people 
highly annoyed versus the outdoor A- weighted day night level. Figure 12 illustrates this dose-response 
relationship for the results from this study. A subjective response rating greater than seven was considered 
a highly annoyed response. The circles in the plot represent the results of this study using the measured 
indoor levels. The shaded area represents the range of data when the indoor levels are transformed to 
outdoor levels by the addition of a 15 to 20 dB correction for house attenuation (ref. 3). Comparison of 
the sonic boom data with the Federal Aviation Administration (FAA) guideline (Shultz curve, ref. 4) 
indicates that a greater percentage of the test subjects were highly annoyed by the simulated sonic booms 
than would be expected for aircraft flyover noise at a given level. 



Figure 12. 
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PERCENTAGE OF TEST SUBJECTS HIGHLY ANNOYED AS A FUNCTION OF C-WEIGHTED 

DAY-NIGHT LEVEL 


The impact of impulse noise has usually been examined in terms of the percent of people highly 
annoyed versus the outdoor C- weighted day night level. Figure 13 illustrates this dose-response 
relationship for the results from this study. A subjective response rating greater than seven was considered 
a highly annoyed response. The circles in the plot represent the results of this study using the measured 
indoor levels. The shaded area represents the range of data when the indoor levels are transformed to 
outdoor levels by the addition of a 15 to 20 dB correction for house attenuation (ref. 3). Comparison of 
the sonic boom data with the Committee on Hearing, Bioacoustics, and Biomechanics (CHABA) 
recommended curve (ref. 5) indicates that a smaller percentage of the test subjects were highly annoyed by 
the simulated sonic booms than would be expected by other impulse noises at a given level. 



Figure 13. 
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CONCLUSIONS 


Figure 14 summarizes the conclusions reached so far from the in-home study of subjective response 
to simulated sonic booms. 


. Proven, reliable system for in-home testing 

. In-home testing requires large resource commitment 

. "Level + 10 * log(n)" confirmed for multiple 
occurrences 

. Perceived level was best annoyance predictor 
. Startle increases annoyance 

Figure 14. 
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INTRODUCTION 


The presentation this morning is about residents’ reactions to sonic booms in a long-term 
sonic boom exposure environment. Although two phases of the data collection have been completed, 
the analysis of the data has only begun. The results are thus preliminary. The list of four authors 
reflects the complex multi-disciplinary character of any field study such as this one. Carey Moulton 
is responsible for Wyle Laboratories’ acoustical data collection effort. Robert Baumgartner and Jeff 
Thomas of HBRS, a social science research firm, are responsible for the social survey field work and 
data processing. The study is supported by the NASA Langley Research Center. 

The study has several objectives. The preliminary data we will consider today address two of 
the primary objectives. The first objective is to describe the reactions to sonic booms of people who 
are living where sonic booms are a routine, recurring feature of the acoustical environment. The 
second objective is to compare these residents’ reactions to the reactions of residents who hear 
conventional aircraft noise around airports. 

Here is an overview of the presentation. This study will first be placed in the context of 
previous community survey research on sonic booms. Next the noise measurement program will be 
briefly described and you will hear part of a social survey interview. Finally data will be presented 
on the residents’ reactions and these reactions will be compared with reactions to conventional 
aircraft. 

Twelve community studies of residents’ reactions to sonic booms were conducted in the 
United States and Europe in the 1960’s and early 1970’s. None of the 12 studies combined three 
essential ingredients that are found in the present study. Residents’ long-term responses are related to 
a measured noise environment. Sonic booms are a permanent feature of the residential environment. 
The respondents’ do not live on a military base. The present study is important because it provides 
the first dose/response relationship for sonic booms that could be expected to apply to residents in 
civilian residential areas. 
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DATA COLLECTION 


The present study was conducted in a Western state in five small communities near military 
operations areas. The communities are separated by from 18 to 104 miles and are from 50 to 150 
miles from the Air Force base from which the training operations are conducted. Less than three 
percent of the residents have employment related to the Air Force base. Although this is a sparsely 
settled area, all of the respondents live in these five settlements and fewer than five percent live on a 
ranch or farm. 

Two phases of the study have been completed. The first-phase interviews were conducted in 
April and May of 1993 and the second phase in December of 1993. In both phases residents were 
asked about sonic booms during the "last six months". Thus, before interviews could be conducted, 
the preceding months’ booms had to be measured. Unattended noise measurements were made with 
Boom Event Analyzer Recorders (BEARs) (Lee, et al., 1989). The BEAR is a 16-bit 
microprocessor-based instrument equipped with a special pressure transducer. The BEAR 
continuously samples the background noise then captures and stores the wave form of loud impulsive 
sounds along with other identifying information. For these measurements the BEARs were set to 
store noise events which exceeded 107 dB, lasted at least 15 msec., had a positive pulse time of 10 
msec., and a rise time of at least 6 dB/35msec. measured just before the peak. The stored events 
were later downloaded and examined to eliminate thunder and occasional other events that did not 
have acoustical profiles that are characteristic of sonic booms. The BEARs were located in 
weatherproof boxes on government property or in a cooperative resident’s backyard at one or two 
locations in each community. 

Depending on the site, approximately 20 to 80 percent of the days in the six months 
preceding the interviews were monitored with functioning BEARs. Late BEAR placement at two 
sites for Phase I and sporadic equipment malfunctions during both phases at all sites resulted in 20 to 
80 percent monitoring rates rather than 100 percent monitoring rates. These data were used to 
calculate average daily exposures for each of the five sites during each of the two 6-month periods. 
These exposures vary from an average of one boom per day to one boom per 10 days and from a 
Day-night Average Noise Level (DNL) of 41 dB(A) [58 dB(C)] to 20 dB(A) [36 dB(C)]. 
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VERY MUCH 
MODERATELY 
A LITTLE 
NOT AT ALL 


Figure 1 

: Response card for noise questions 


Q6. We want to learn how you feel about the neighborhood right around here and 
about any advantages that make you feel it is a good place to live* In the 
six months since we last talked to you what are the one or two things you 
have liked most about this area? 


Q7. How about any things that are disadvantages. What are the one or two 

disadvantages that you have disliked the most about this area in the last 
six months? 


Figure 2 % Beginning page of sonic boom questionnaire 

The residents’ reactions were obtained through fixed-format, interviewer-administered, face- 
to-face questionnaires. The questionnaire was developed during three rounds of pretests at three 
other locations in the United States at which booms are heard. The pretests provided the basis for 
improvements such as reducing unneeded material (all respondents were familiar with the term 
"sonic booms from jets") and maintaining rapport despite detailed questions about sonic boom 
effects. Interviewing procedures and the questionnaire were designed so that respondents were not 
initially alerted to the subject matter of the survey. 

The professional interviewers on this project received an additional day of study-specific 
training for this study. Question-specific instructions were prepared. Control was maintained over 
interviewing procedures through on-site supervision, interview tape recordings and verification of 20 
percent of the interviews by staff visits or follow-up telephone calls. An overall response rate of 78 
percent was achieved. 

Homes were randomly selected in three communities with strict procedures for identifying 
pre-selected respondents within homes. In the two smallest communities all eligible adults were 
interviewed. Due to the small study populations, about 400 of the Phase I respondents were 
reinterviewed during Phase II with a slightly shortened questionnaire. 
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Now some questions about noises you might have heard when you have been at home in 
the last six months. 

a. What are some of the different types of noises you have heard around here? 

(PROBE: Anything else?) [MARK "VOL" FOR VOLUNTEERED NOISES] 

b. [ASK FOR ALL NOISES NOT VOLUNTEERED] In the last six months, have you ever 
heard the noise from ...(cars or trucks or other road traffic going by)... 
when you were here at home? 

!!: COMPLETE ENTIRE LIST WITH b BEFORE STARTING C] 

c. Here is an "AMOUNT" card for choosing your answer for the next question. [HAND 
CARD A TO RESPONDENT] [ASK FOR EACH SOUND HEARD] 

During the last six months has the noise from ...(cars or trucks or other road 
traffic going by ).. .bothered or annoyed you very much, moderately, a little, 
or not at all? 



c. BOTHERED OR ANNOYED 


VERY 

MODER- 

A 

NOT AT 

DK 

MUCH 

ATELY 

LITTLE 

ALL 


(i) 

(2) 

(3) 

(4) 

(8) 

VERY 

MODER 

LITTLE 

NOT 

DK 

1 

2 

3 

4 

8 


Cars or trucks or other road 
traffic going by 


ii. Motorcycles 


iii. Neighbors' tools or outdoor 
equipment 


iv. [REPEAT FULL QUESTION] 
..Sonic booms from jets 


Any other explosions or 
bangs or booms (besides the 
sonic booms) (DESCRIBE) 


Low-flying jet aircraft 


vii. Any other airplanes 

(besides the low-flying 
jets) 

(DESCRIBE) 


viii. [DESCRIBE ANY OTHER 
VOLUNTEERED NOISES HERE] 


1 VOL.JI 

2 YES JL 



1 VOL . 

2 YES | 


3 NO } 

4 DK | 


1 VOL. 

2 YES | 


3N0 | 

4 DK jj. 


1 VOL.U 

2 YES | 


3 NO | 

4 DK | 


1 VOL.JJ 


NOT 

DK 

MARK X 

4 

8 

AT Q12 

[MARK 

XS =>] 


& Q20 


VERY 

MODER 

LITTLE 

NOT 

DK 

1 

2 

3 

4 

8 



VERY 

MODER 

LITTLE 

NOT 

DK 

1 

2 

3 

4 

8 



Figure 3: Second page of sonic boom questionnaire 

I have a tape recording of one of the interviews for you to listen to now. The interviewer 


212 





























































switched on the recorder after entering the home, making the routine selection of the respondent, and 
obtaining permission to turn on the tape recorder. The interview begins with Question Q6. (Items 
Q1 to Q5 were for sample identification.) For the respondent this is almost entirely an aural event. 
The respondent sees only the interviewer and the occasional answer card (Figure 1). Only the 
interviewer sees the questionnaire, the relevant parts of which are reproduced in Figure 2 and 
Figure 3. 

The sonic boom question (Question 8.iv in Figure 3) exactly parallels the aircraft noise 
question that has been previously used with over 10,000 residents near airports. The remainder of 
the questionnaire obtained reactions on other sonic boom annoyance scales and obtained detailed 
information about vibration, startle reactions, suspected damage, types of activity interference, 
demographic characteristics, and more general attitudes. 


RESULTS 

After the 1042 interviews were entered into a computer file and merged with the noise data, 
the reactions in each community were plotted by noise level. Figure 4 shows that at, for example, 
about DNL 25 to 35 dB(A), about 30 percent of the respondents said that they were "very much" 
annoyed by sonic booms (Question 8.iv. in Figure 3). To obtain a broader perspective on the 
meaning of that degree of annoyance, the respondents chose numbers to compare their feelings about 
sonic booms (arbitrarily set at a score of 100) with their feelings about other life-familiar events for 
which they could choose any number. Logarithmic averages of the ratios of their numerical ratings 
of sonic booms and other events indicate that the 79 percent of the respondents in the 25 to 35 dB 
range who expressed any annoyance with sonic booms felt that sonic booms were about as annoying 
as "hearing big noisy trucks if you lived near an intersection" or "having a dog next door that 
regularly barks in the middle of the night". As for all noise sources there is considerable variation 
in the respondents’ reactions. In the 25 to 35 dB range, for example, about 21 percent are "not at all 
annoyed" by sonic booms. The remaining 79 percent who are at least a little annoyed consist of 27 
percent who are "a little annoyed", 22 percent who are "moderately" annoyed, and 30 percent who 
choose the highest category of "very much" annoyed. 

From a regulatory perspective the most important comparison is not with reactions to 
"barking dogs" but with reactions to other aircraft. In later analyses it is planned to relate the sonic 
boom results directly to 8 surveys that used the annoyance scales that are included in the sonic boom 
survey. The present preliminary results can be directly compared with one of those surveys; a 1985 
survey of residents around five airports in the United Kingdom. Figure 4 presents this comparison 
for DNL (A). The figure indicates that the level of annoyance that is expressed in the sonic boom 
study at about 40 dB is not reached in this conventional aircraft noise study until about 60 to 70 dB. 

The next figure (Figure 5) adds another line. The new line is the dose/response relationship 
predicted by a logistic regression analysis of 400 data points from 26 community noise survey data 
sets (Finegold, Harris, and von Gierke, 1994; Federal.., 1992). These 400 data points were drawn 
from a larger set of 453 data points that appeared in a dose/response synthesis study (Fidell, Barber, 
and Schultz: 1991). The new line is thus broadly consistent with the previous conclusion that sonic 
booms are much more annoying than other community noises. 
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Figure 4: Comparison of Western Boom Study and UK 5-Airport Study 
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Figure 5: Comparison of Boom, UK 5-Airport and FICON guideline 


In the following figure (Figure 6) the sonic boom sites are compared to all 453 data points 
from the complete synthesis study (Fidell, Barber, and Schultz: 1991). The general conclusion 
remains unchanged; sonic booms appear to be more annoying than would be predicted by other 
community noise data. While all of the displays of data come to this same general conclusion, 
several difficulties interfere with drawing more precise conclusions about the size of the difference 
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between reactions to sonic booms and to other aircraft. First, the 453 data points are based on a 
diverse set of annoyance questions, only one of which was the same as was used in the sonic boom 
survey. Second, the 453 data points represent studies of reactions to a diverse set of noise sources 
including aircraft, road traffic and undifferentiated community noise. Third, the 453 data points 
from the 29 studies differ from one another in more fundamental ways that hinder any comparison. 
The points marked with open symbols in Figure 6, for example, are based on data from studies with 
the following types of weaknesses: 1) the respondent was not asked about the measured noise 
environment, 2) the annoyance measure is relative to local features that vary between sites, (3) a 
transportation noise source is not specified in the question, or 4) the annoyance measure is a complex 
combination of questions that is not clearly definable as high annoyance (Fields, 1994). 


Figure 6 raises one other issue; there is very little overlap between the noise levels observed 
for the sonic boom study and the noise levels observed for the other surveys. Although the three 
lowest overlapping data points in the synthesis study (between DNL 29 and 43 dB) used the same 
annoyance scale as was used in the boom survey, all three come from a railway noise survey in 
which the respondents have been shown to be less annoyed than respondents to conventional aircraft 
noise surveys (Fields and Walker, 1982). These very low sonic boom exposures are a result of the 
very low frequency of the sonic boom events. The use of energy averaging for such infrequent 
events has not been examined in previous studies. If it were speculated that respondents at such low 
levels do not react to the long-term average, but rather to the worst sonic boom days, then the sonic 
boom reactions should be plotted at higher noise levels. To examine the implications of such 
speculations, the sonic boom reactions in Figure 7 are plotted at noise levels which are DNL 14 
dB(A) higher than previously. The 14 dB difference is approximately the average difference between 
the average noise level for all days and for the three highest days. Even after this speculative 
adjustment that does not make a corresponding adjustment for the three highest days around 
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Figure 7: Speculative comparison if respondents misunderstood and reported reactions to the 


worst three days for only the sonic boom study 
conventional airports, the sonic boom reactions are stronger than those to other noise sources. 



The final figure (Figure 8) shows four logistic regression lines; one for each of the two boom 
study phases, one for the five-airport UK study and one for the 400 point synthesis. This figure 
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raises at least one additional point: the difference between the two phases of the boom study. While 
we can presently theorize about the reasons for such a difference, the data have still to be examined. 
The precision of the noise estimates and the dose/response relationships must both be examined 
before these and other issues can be addressed. 

It is important to conclude by emphasizing that these are preliminary results from an on-going 
project. An additional data collection phase is planned, some aspects of the noise data have still to 
be examined, and detailed analyses of the social survey data have not begun. The present data 
indicate that sonic booms are more annoying than would be predicted from the findings from one 
carefully matched aircraft noise survey and from widely-accepted summaries of dose/response 
relationships. More analysis is planned before the size of this difference is quantified and before 
explanations for the differences in reactions can be examined with these data. 
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